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THE  MISSION  OF  ACARD 
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PREFACE 


Ablating  heat  shields  have  successfully  protected  many  hypersonic  vehicles  during 
atmospheric  re-entry.  Detailed  knowledge  of  the  characteristics  of  many  materials  and  of 
the  interaction  mechanisms  between  materials  and  high-temperature  environments  exists. 

Using  this  information,  very  reliable  passive  heat  shields  can  be  designed. 

There  is  an  extensive  literature  containing  the  results  of  both  analytical  and  experi¬ 
mental  investigations  of  ablation.  Although  the  technology  has  reached  a  mature  status, 
new  materials  arc  continually  being  developed  and  passive  ablating  systems  are  of  great 
interest  to  heat-shield  designers.  The  well-developed  literature  and  continuing  interest 
in  ablation  technology  make  this  a  particularly  suitable  subject  for  an  AGARDograph.  The 
intimate  coupling  between  the  fluid  mechanics  of  high-temperature  boundary  layers  and 
ablating  materials  makes  sponsorship  by  the  Fluid  Dynamics  Panel  very  appropriate. 

This  AGARDograph  review  the  advances  in  ablation  technology,  the  present  status, 
and  indicates  possible  future  directions.  Ablation  analyses  are  discussed  starting  with  the 
original  simple  model  in  which  the  amount  of  material  removed  is  computed  from  the  heat 
transfer  to  a  non-ablating  body  and  the  effective  heat  of  ablation.  Analytical  methods  are 
carried  through  the  development  of  the  most  sophisticated  approach  which  considers  the 
coupling  between  the  boundary  layer,  the  ablating  materials,  and  includes  complete  non¬ 
equilibrium  chemistry.  The  behavior  of  various  materials  is  compared.  The  status  of  today’s 
experimental  facilities  used  to  determine  ablation  behavior  is  reviewed,  including  test  tech¬ 
niques  and  instrumentation.  Finally,  the  shape  changes  and  changes  in  surface  roughness 
which  result  from  ablation  are  discussed.  These  changes  obviously  interact  with  vehicle 
performance,  flow  fields,  and  aerodynamic  characteristics. 

The  authors  have  covered,  with  remarkable  completeness,  a  large  complex  subject.  The 
bibliography,  as  well  as  the  material  covered  here,  should  be  of  considerable  use  to  those 
interested  in  ablation  technology  and  the  design  of  ablating  heat  shields. 

Dr.R.E.  WILSON 

Fluid  Dynamics  Panel  Member 


The  AGARD  Fluid  Dynamics  Panel  wishes  to  thank  Dr.  Wilson  for  his  contribution  in  under¬ 
taking  the  review  and  editing  of  this  AGARDograph. 
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SUMMARY i 


Two  decades  of  advances  in  ablation  technology  are  reviewed  vlth  particular  emphasis  on 
developments  of  the  last  eight  years.  Leveling  off  of  research  activities  is  noted  reflecting  both  the 
apace  systems  requirements  and  acceptance  of  ablation  as  an  established  meant  of  thenoal  protection. 

Physl co-mathematlcal  approaches  to  analysis  of  ablation  processes  are  brought  out  rsnglng 
from  simplistic  heat  of  ablation  concept  to  the  sophisticated  aialysia  of  laminar  and  turbulent  flow,  non 
equilibrium  ablation  In  chemically  reacting  flow  fields  Including  surface  reactions  and  radiative  heat 
transfer.  Particular  emphasis  Is  given  to  the  eynargintlc  mechanisms  occurring  during  the  ablation 
process.  Examples  of  nosetlp  interaction  with  eero/thermodynsmic  environment  at  the  surface  and  in  depth 
phenomena  la  given,  as  veil  a e  strlation  phenomena  on  the  afterbody  ere  discussed. 

Graphite  Is  used  as  an  example  of  synergistic  ablation  phenomena  and  effects  because  of  the 
interest  noted  In  the  literature  and  because  o'.ner  materials  may  be  treated  In  a  similar  fashion. 

Graphite  eroalon  mechanisms  and  Its  microstructure  characterisation  are  reviewed.  The  effects  of  fillers 
binders,  porosity,  and  other  parameters  are  discussed  as  well  as  the  synergistic  effects  of  the  above 
paraswters  are  also  approached  In  a  novel  fashion. 


Material  property  and  environmental  effects  on  material  performance  ere  also  reviewed  as  well 
as  the  test  simulation  requirements.  Requirements  for  further  etudiee  of  mechanical  erosion,  other 
synergistic  affects,  and  extension  of  facilities  range  of  performance  are  also  noted. 


LIST  OF  SYMBOLS 

s 


C 

Pf 


c 

pl 


etmoe 


P5 


Local  (dimensionless)  pure  convective  heat  tranefer  coefficient  or 
Stanton  numbers  including  dissociation,  mass  Injection  and  surface 
reactions 

Lccal  (dimensionless)  diffusion  (mass  transfer)  Stanton  number 


Local  (dimensionless)  cold  wall  heat  transfer  coefficient  or  Stanton 
number  uncoupled  freo  species  solution 

Specie  mass  fraction  or  concentration  of  ith  specie,  Pt/P 

Specific  heat  for  composite  mixture  of  boundary  layer  species  at 
constant  pressure 

Specific  heat  for  1th  species  at  constant  pressure 


Specific  hast  of  condensed  phase  in  the  decomposing  sbletor 


Specific  heat  of  mixture  of  bourdary  layer  species, 

Specific  heat  for  composite  mixture  of  gaseous  products  of 
depolymerization  at  constant  pressure 
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Activation  energy  uaed  in  Arrhenius'  rate  expression  for  the 
representative  solid  phaae  reain  component 

Activation  energy  used  in  Arrhenius'  rate  expression  for  the 
representative  completely  volatile  resin  component  o12 

Activation  energy  used  in  Arrhenius'  rate  expression  for  the 
representative  completely  volatile  nylon  reinforcement  o? 

Activation  energy  for  Arrhenius '  rate  expression  describing  rate 
limited  surface  erosion  or  removal  by  oxidation  of  graphite  Cj 

Activation  energy  used  in  Arrhenius'  rate  expression  describing  a 
form  of  deposition  interaction  or  coking  for  the  ablation  model 

Ratio  of  char  density  to  virgin  plastic  density  of  model  ablator 

Number  fraction  of  particles  eroded 

Maes  fraction  of  particles 

Solid  fraction  mechanically  eroded 

Melt  or  evaporation  fraction,  m^/m^,  giving  the  ratio  of  the  rate  of 
evaporation  to  the  total  ablation  rate  in  the  liquid  layer  equation 

External  radiation  potential  from  gloving"  internal  ablation 
reactions  as  a  result  of  char  transparency 


Reduced  (dimensionless)  total  enthalpy  function,  refers  'o  I/I( 
in  the  laminar  boundary  layer  solutions 

Reduced  (dimensionless)  total  "sensible"  er.thalpy  function, 

I  /(I  )  in  the  laminar  boundary  layer  solutions 
f  f  e 

r  ,  f  ^  o » 

Enthalpy  function  defined  m  /  ^  (  Cjdt  ♦  / 

1  *  o 


/  A 

Sensible  enthalpy  function  defined  »«  [  C.  I  C  dt 

1  ‘o  Hi 

fT  o 

Specie  enthalpy  function,  I  C  dt  e  h^ 

■lo  P1 

Chemical  enthalpy  of  formation  of  the  ith  specie* 

Combined  heat  of  melting  «nd  vaporltation  (also  include*  eombu»tlon 
and  internal  reaction  effect*  if  not  othervise  accounted  for)  in  the 
liquid  layer  expression 

Wall  enthalpy  (at  ablator-boundary  layer  interface) 

Lumped  enthalpy  of  ma»s  transfer  and  vaporitation  or  sublimation 
term 

Total  enthalpy  at  boundary  layer  edge ,  same  »a  lg 

f  t  o  2 

Total  enthalpy  function  defined  as  J  ( I  C  dt  ♦  ) ♦  u  /? 

i  *  ** 

'otal  "sensible"  enthalpy  function  defined  a*  [  Cj  Cp^  t  u?/£ 
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Variable  Id  the  "gloving"  expression 

Frequsney  factor  used  In  Arrhenius'  rat*  sxprssslon 
Also  exponent  of  y  to  signify  tvo  dimensional  (K*c)  or  axi- 
symmetrlcal  (k»l)  6 
Uao  volume  fraction 

Frequency  coefficient  ua«*d  is  Arrhenius’  rat*  expression  for  the 
representative  tolld  phase  resin  component  0^. 

Frequency  coefficient  uaad  in  Airhenies'  rat*  expression  for  th* 
representative  completely  volatile  reein  component  o, 

Frequency  coefficient  ueed  in  Arrhenlue'  rat*  expression  for  the 
representative  coaplstely  volatile  nylon  relnforceseot  pj 

Frequency  coefficient  used  in  Arrhenius'  rate  expression  describing 
a  fora  of  deposition  interaction  or  coking  for  the  ablation  model 

Max*  fraction  of  oxygen  at  th*  outer  edge  of  the  boundary  layer 


Binder  mass  fraction 


Geometric  mean  nast  flux 

Mass  flux  due  to  surface  erosion  or  removal  by  oxidation  of  graphite 
C  when  in  the  rate  limited  temperature  regime 

Maaa  flux  due  to  solid  surface  erosion 

Composite  maaa  flux  of  gasecua  products  of  decomposition 

Mass  flux  of  ablation  gas  products.  Includes  both  m^,  and  n^„ 

Kaee  flux  of  solid  species  of  surface  material 
Mass  flux  of  lndepth  pyrolysis  gas  species 
Total  injected  rasa  flux  fr„s,  the  ablating  surface 
Mass  flux  due  to  mechanical  erosion  (liquid  or  solid) 

Mean  molecular  veight  of  boundary  lays-  species 


Mean  molecular  .eight  of  injected  species 
Molecular  veight  of  atomic  oxygen 
Molecular  veight  of  molecular  o>ygen 

rrceiure  exponent  appearing  In  rate  lav  expression 

Reaction  order  exponent  used  In  rate  lav  expression  for  the 
representative  solid  phase  resin  component 

Reaction  order  exponent  used  In  rste  lav  expression  for  the 
representative  completely  volatile  resin  component  o^g 


it 


n,  Reaction  orler  exponent  used  in  rete  lev  expression  for  the 

representative  completely  volet lie  pylon  reinforcement  pg. 

Alio  reectlon  order  exponent  for  oxldetlon  of  graphite  by  0  In 
rete  lev  expreieion 

n,  Reectlon  order  exponent  for  oxldetlon  of  graphite  by  0.  In  rrte 

3  lev  expreieion 

n  Local  partial  preeeure  of  char  depoeltion  exponent  In  coking  rete 

c  lev  expreieion 

n^  Reectlon  order  exponent  for  a  rate  lev  expreieion  describing  a  font 

of  deposition  for  the  ablation  model 

p  Local  pressure 

pc  Local  partial  preeiure  of  coking  char  materiel 

Local  boundary  layer  edge  preeiure 

p*  Oenerelixed  preeiure 

P  Molecular  Prandtl  Number 


Eddy  (turbulent)  Prandtl  Number 


^blocking 


beat  transfer  rete  contribution  due  to  "heat  blockage"  by  mass 
transfer 


^chemical 

s> 


*r 


S-ad 


Heat  transfer  rate  contribution  due  to  chemical  reactions  in 
surfece  energy  balance  equation 

Cold  vail  pure  convective  heat  transfer  rate  evaluated  by  Eckert's 
reference  enthalpy  method. 

Also  stagnation  point  heat  transfer  rate 

Radiative  plus  convective  heat  transfer  rate 


Net  radiative  heat  transfer  rete  into  surface  q  -q  or 


-t  of  ♦  a  q 
v  v  v  Tad 


Radiative  heat  flux  to  surface  from  environment 


q 


j 


Net  heat  transfer  rate  or  flux  through  surface  to  the  Interior  of 
ablator  *  -(x  dt/dy)^ 


<v 


r 

R 


R 


o 


R 


y 


Aerodynamic  heat  transfer  rate  to  the  surface  of  the  ablator  from 
the  boundary  layer 


1/2 

Recovery  factor,  *  P  (Tor  lam.  B.L.),  * 
Universal  gas  constant 

Reflectivity  of  the  char  at  the  surface  of 
"gloving1'  term 

Reflectivity  of  the  char  material  in  depth 
evaluation  of  df/dy  appearing  In  "gloving" 


P1/3  (for  turb.  B-L.) 


the  ablator  appearing 


st  location  of 
term 


In 


e 


Steady  state  recession  rate 


t 


Time 
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Binder  bridge  thicksets 
Local  temperature 

Local  "vail"  or  surface  temperature 

Streamvite  component  of  flov  field  velocity  vector  v  parallel  to  s 
direction 

Transverse  component  of  flov  field  velocity  vector  v  normal  to  s 
direction  and  parallel  to  y  direction 

Annihilation  or  creation  of  species  (i.e.,  mass  rate  of  change  of 
specie)  through  chemical  reactions  vithin  the  boundary  layer 

Special  curvilinear  coordinate  normal  to  s  coordinate  and  transverse 
to  ablator  surface,  also  depth  into  material 

Transformed  depth  into  miterial  (Eq  '*!) 

Absorptance  in  depth  used  in  "gloving"  expression. 

Also  viscous  resistance  coefficient 


Absorptance  of  surface  of  the  ablator 


Laminar  and  turbulent  boundary  layer  mess  injection  parameter 

Specific  heat  ratio,  CWC^,  Tor  ideal  gaa;  also  given  as  the  char 
mass  fraction  of  pyrolyiing  material 

Volume  ‘‘ractlon  of  a  degradation  polymer  in  a  virgin  ablr.tor 
composite  (represented  as  Phenolic  Nylon  for  the  ablator  model) 

Char  depth  measured  from  the  outer  surface  of  an  ablation  test 
model 

Ablator  nodal  tblcKneas  used  in  transient  heat  conduction  solution 

Enerror  released  by  combustion  reactions  between  the  gaseous  surface 
material,  pyrolysis  products  and  environment  combined  Into  one  term 

Boundary  layer  enthalpy  potential  for  the  free  stream  gas 

Enthalpy  of  deposition 

Enthalpy  of  depoiymerisatlon  and  pyrolysis 


Energy  absorbed  by  pyrolysis 


Net  energy  of  decomposition  plus  sensible  enthalpy  vith  original 
temperature  and  state  as  reference 

Porosity 

Ablator  surface  emlsslvlty 
Parameter  defined  by  Eq  36 

Similarity  transformed  (Lee’*  Dorodnitsyn)  or  reduced  body 
coordinate  y  used  in  laminar  boundary  layer  solutions,  also 
transpiration  coefficient  used  in  liquid  layer  equation 

Viscosity 

Molecular  (for  gas  specie*)  thermal  conductivity 
iddy  (turbulent)  thermal  conductivity 


Parameter  defined  by  Eq  39 

Local  composite  density  of  all  species  (gaseous-both  boundary  layer 
and  pyrolysis  products) 

Instantaneous  condensed  phase  density  of  the  representative  solid 
phass  resin  component  of  the  model  ablator 


I 


6 


o 


12 


dP 


(pv) 

rp 

(pv)w 


Instantaneous  condensed  phase  density  of  the  representative 
depolymcrlted  end  completely  volatile  recln  coupon -nt  of  the  model 
ablator 

Instantaneous  eondenaed  phaac  density  of  the  repraaentatlve 
depolyaerlted  and  completely  volatile  nylon  reinforcement  component 
of  the  model  ablator 

Inatantaneoua  eondenaed  phaae  density  due  to  a  form  of  deposition 


Instantaneous  velghted  sum  of  local  polymer  condensed  phaae 
densities  (o  Dj) 

Density  of  the  1th  species  In  the  boundary  layer  of  the  same  unit 
volume  as  p 

Density  of  original  ablative  material  before  decomposition 


Total  Instantaneous  condensed  phase  densities  of  local  polymer  and 
deposition  contributions 

Hass  flux  of  pyrolysis  gsses 


Total  ablation  mass  flux 


o 


Stefan-Boltzmann  constant 


o' 


Eroded  surface  depth 


T 


v 


Aerodynamic  shear  stress 


Subscripts 

b 

e 


f 

1 

k' 


k" 


J.  k,  1 

o 

r 

r 

P 

s 

th 

v 

V 


Binder 

Denotes  boundary  layer  edge  value. 

Also  end  effect  region  properties 

Denotes  sensible  enthalpy  term 

Denotes  ith  species,  initial  value;  or  1th  nodal  laminar 

Denotes  ablator  char  solid  phase  species 

Denotes  ablator  pyrolysis  ks»  species 

Refers  to  orthogonal  rectangular  coordinates  (3  dim.) 

Refers  to  stagnation  point  value  or  virgin  material 

Denoter  recovery  value 

Pyrolysis  product 

Denotes  condensed  phase;  or  net  surface  value 

Thermochemical 

Pore  volume 

Denotes  "vail"  or  ablator  surface  value 


Denotes  free  stream  value. 
Also  thick  sample  properties 
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I. 0  UTRODUCTION 

X*  the  aeientiflc/engineering  community  applauded  the  success  of  the  early  unmanned  and  manned 
reentry  flight!,  it  recognlied  the  (mergence  of  ablation  ti  a  reliable  neant  of  therual  protection.  Theae 
successes  lntanalfied  the  reaeerch  efforta,  vfcich  over  a  decade  and  a  half  have  advanced  the  analytical 
toola  froa  the  relatively  simple  "beat  of  ablation"  approach  (1)  to  aophiatieated  computer  codeo  which 
allov  e  detailed  description  of  toe  ablation  proeeaaea  (2, 1,1,5).  Parallel  with  the  analytical  develop- 
aenta  were  the  material  formulation*  required  to  tatitfy  the  more  specialized  and  usually  more  stringent 
mission  requirements.  These  combinations  produced  a  new  challenge,  it  was  no  longer  sufficient  to  be 
able  to  predict  the  gross  behevlor  of  ablation  materials.  It  uaa  nov  required  to  understand  the  oicro- 
ecoplc  behavior  in  order  to  adequately  give  direction  to  the  nev  materials  development  required  to  satisfy 
the  ever  increasing  system  demand,  furthermore,  the  efficiency  of  the  very  process  of  ablation  (removal 
of  mass  by  theraochemlcal  or  thermomechanical  proceeaet)  had  to  ti  maximized  and  in  the  limit  reduced  to 
perhaps  just  the  surface  oxidation  at  in  the  case  of  graphites  in  some  environments. 

The  purpose  of  this  AGARDograph  is  to  review  the  advances  lb-  and  the  statua  of-  the  ablation 
technclogy  after  nearly  two  decades  of  lta  application  to  various  space  reentry  and  propulsive  system*. 

Many  review  papers  (with  extensive  bibliographies )  heve  been  published  during  this  period  (6,  7,  8,  9,  10, 

II,  12)  each  with  somewhat  different  emphasis  on  various  aspects  of  ablative  phenomena.  Some  of  the 
reviews  stressed  ebl-  tion  mechanisms  or  materiel  behevlor,  while  others  treated  environment  Interaction 
with  ablative  system  end  experimental  simulation  of  ablation  phenomena  in  flight.  In  the  latter  part  of 
tha  alxtlas  production  of  revlev  papara  nearly  matched  the  publication  of  papers  treating  new  aspect*  of 
ablation  phenomena.  This  indicates  a  leveling  off  of  research  activities  which  may  be  due  to  aatls factory 
solution  of  engineering  problems  associated  with  present  or  presently  contemplated  space  missions  -  or  in 
other  words  attainment  of  maturity  by  ablation  which  now  has  Income  an  established  r»ans  of  thermal  pro¬ 
tection  and  found  its  way  into  handbooks.  This  need  not  necessarily  mean  a  halt  in  research  end  develop¬ 
ment  activities,  but  it  does  orlant  such  activities  towards  specific  mission  requirement*  rather  then 
broad  theoretical  studies  so  typical  of  the  fifties  and  early  sixties.  This  decreased  emphasis  cn  ablation 
research  per  as  may  be  caused  by  the  fact  that  advanced  missions  In  many  cases  operate  either  la 
environments  which  invite  etuis,  of  nan-receding  protection  systems  or  demand  reuse  of  the  surface  with 
minimum  of  refurbishment.  The  attractiveness  of  graphites  or  low  density  polymeries  cannot  be  underesti¬ 
mated  ard  thus  studies  of  "near”  ablation  phenomena  take  precedence  over  ablation.  Still,  ablation  is 

on*  of  the  swat  affective  and  economic  means  of  protection  of  reentry  vehicles.  X  previous  AOARbcgraph  (8) 
revlev  of  the  related  subject  of  aero/thermoebemistry  of  ablation  by  one  of  the  authors  of  this  >ubllca- 
tloo  carried  the  advances  is  ablation  into  the  sixties  and  will  be  treated  here  as  e  departing  ;  Dint  to 
avoid  rspatitlon.  Also  the  authors  have  contributed  two  chapter*  (13,  1U )  on  the  subject*  of  S'  Lation 
and  thermal  protection  systems,  respectively,  which  will  be  used  here  where  required,  frlven  b:  require¬ 
ment*  of  various  space  missions  briefly  noted  above,  the  developments,  scant  that  they  have  bee  in  the 
laat  five  years,  concentrated  in  the  following  areas: 

1.  Internal  pyrolysis  phenomena  for  lov  density  ablator*. 

2.  Behavior  of  graphites  in  high  pressure  and  temperature  environments  including  internal 
pyrolysis  and  mechanical  erosion, 

3.  Fmplrical  evaluation  of  mechanical  erosion  of  ablators  and  analysis  of  shape  c.  angea 
due  to  thermochemical/thermomechanical  forces. 

U,  Response  of  materials  in  environment  with  high  radiative  fluxes,  and  including  presence 
of  massive  ablation. 

5.  Better  understanding  of  surface  and  boundary  layer  chemistry. 

6.  Analysis  of  various  synergistic  mechanisms  Including  attempts  on  treatment  of  coupled 
lasctlng  boundary  layers  with  radiative  heat  tranefer,  non-equilibrium  ablation,  and 
surface  end  flow  rield  chemistry  for  laminar  and  turbulent  flow. 

7.  Design  and  operation  of  improved  are  facilities  capable  of  very  high  pressures  (-200  atm) 
and/or  radiative  heat  fluxes  in  combination  with  previously  attained  high  enthalpy 

plasmas . 

The  dlacusalon  of  the  environment  per  se  la  out  of  the  scope  of  this  paper.  The  reader  it 
referred  to  (15,  l6)  for  euch  review.  However,  the  interaction  of  environment  with  material  response  and 
vehicle  performance  is  critical  in  evaluation  of  ablation  phenomena  and  as  such  will  be  reviewed  here. 

It  1*  hoped  that  this  AGARDograpb  with  it*  review  of  the  etate-of-the-art  in  ablation  analysis,  materials 
behavior  tasting,  and  study  of  synergistic  mechanisms  will  be  of  help  to  the  scientist  and  engineer 
studying  and  designing  apace  vehicles  and  developing  material*  used  for  thermal  protection. 

Further  work  in  extending  the  understanding  of  synergistic  phenomena,  mechanical  erotica:  and 
extending  the  range  of  experimental  simulation  facilities  is  still  in  order,  as  the  state-of-the-art  us 
still  not  sufficient  for  some  of  the  future  Interplanetary  missions  which  cay  be  contemplated. 

2.0  ABLATION  ANALYSES 

2.1  General  Considerations 

The  environmental  parameter.,  associated  with  hypervelocity  flight  in  planetary  atmospheres 
affecting  ablation  mecbanltm*  ore:  gaa  temperature,  enthalpy,  rode  of  energy  transport,  total  heat  load, 
peak  beating  rate,  shape  of  the  heat  pulse,  duration  of  the  pulse,  external  and  internal  pretsure  forces, 
aerodynamic  (hear,  particle  impingement,  mechanical  and  acoustical  vibration,  inertial  arid  dynamic 
forces,  and  the  chemical  reactivity  (17).  Of  primary  concern  is  the  means  of  accommodating  the  energy 
converted  to  heat;  i.e.,  the  thermal  protection  system.  This  le  accomplished  In  three  ways:  (a)  absorption. 
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(b)  radiation,  and  (e)  reduction  of  the  aarodyuamie  heat  prior  to  It  ranching  the  vehicle  (18). 

The  requirement  to  minimise  hod;  weight  haa  often  led  to  the  choice  of  materials  purposely 
designed  to  undergo  internal  endothermic  reectlane  lending  to  the  format ion  of  a  earhonauceoue  char  when 
heated  and  to  vaporlie  at  the  eurface  exposed  to  the  oypersoalc  environment.  The  energr  absorption 
aatoclated  vlth  eucb  proceieet,  together  vltn  t  *  reduction  of  the  heat  trenafer  coefficient  aa  a  reeult 
of  injection  of  the  producte  of  degradation  and  the  vaporized  epeclee  Into  the  boundary  layer  have  proved 
to  be  extremely  effective  In  providing  reliable  thermal  protection  eyetema  (19)* 

it  haa  been  found  neceeeary  (20)  for  efficient  design  to  have  an  --.curate  mathematical 
description  of  the  complex  heat  and  mass  transfer  phenomena  aatoclated  vlth  such  processes  Inside  the 
material  end  an  efficient  numerical  procedure  for  Implementing  these  calculation*  baaed  on  the  mathematical 
model.  In  practice,  the  ablation  phenomena  are  most  often  encountered  In  situations  vhere  the  environmental 
parameters  are  nonlinear  function*  of  time.  This  fact,  coupled  with  the  heat  shield  physically  being  of 
finite  thickness  makes  it  necessary  to  perform  transient  rather  than  steady  state  calculations  (19).  The 
fully  transient  solution  includes  the  decoeipoaltlon  la  depth  with  the  attendant  density  end  temperature 
gradient*.  Internal  transpiration  and  radiation  gtneratlon,  material  rederoeltlon ,  source  terms  for 
convective  end/or  radiative  Input,  and  provide*  for  a  variety  of  boundary  conditions  or  »urfec*  phenomena 
(20).  With  the  aid  of  present  day  computer*  the  solution  is  practical  even  for  quit*  complex  system*. 

It  haa  been  shown  (19,  21,  22)  that  an  Intricate  relation  between  the  environment  and  >.».t*rl*l 
response  indeed  exists  end  that  the  solution  of  the  problem  must  account  for  them. 

Although  the  transient  problem  can  be  solved  with  a  certain  degree  of  reliability  fur  certain 
combinations  of  environment  and  material*,  or  for  steady  state  experiments,  e  steady  state  solution  Is 
deemed  adequate  and  Is  lee*  expensive  In  terms  of  computing  time.  Often,  in  preliminary  design  phase*,  a 
combination  of  transient  and  steady  state  solutions  is  used  -  causa  the  costs  associated  with  a  comprehen¬ 
sive  transient  analysis  cannot  be  Justified. 

One  euch  condition  where  a  Heady  state  analysis  (i.e.,  neglecting  indepth  phenomena)  can  be 
effectively  utilized  le  In  the  investigation  of  the  ablation  requirement*  for  lifting  vehicles  which  enter 
from  orbit  and  sustain  long  constant  altitude  trajectories  end  minimal  deceleration* .  Here  the  ablative 
process  remain*  in  the  dlffualon-llsilted  regime  and  involves  the  least  Intricate  of  mechanism*.  In  low 
heating  environments  (i.e.,  heat  fluxes  lea*  then  75  Btu/ft2-»ee) ,  the  eurface  recession  has  been  shown 
experimentally  to  be  lnalgniflcant.  The  formation  of  the  cher  and  the  influence  of  pyrolyala  on  material 
properties  may  also  be  considered  Inconsequential  vlth  respect  to  a  material  Insulating  capability.  An 
ae-jumptloj.  of  this  nature  inherently  postulates  an  energy  balance  of  the  contradicting  effects  of  pyrolyal* 
whereby  blockage  or  heat  transfer  by  gaseous  mass  flux  is  compensated  for  by  neglecting  the  Increase  of 
*herral  conductivity  due  to  char  formation.  Hence,  based  on  these  premise*  In  low  heating  environments, 
a  simplified  steady  state  model  nay  he  utilized  for  estimating  ablation  requirements.  Other  conditions 
wnere  similar  approacuea  can  be  taken  to  reduce  the  problem  to  a  steady  state  so)"* ion  an»  dlacuaaed  In 
(20.  23,  2M. 

2.2  Surface  Mcchapiema  -  Transient  Ablation 


The  sevetal  energy  transfer  mechanisms  between  the  boundary  layer  gas  end  the  ablating  surface 
are  shown  In  Figure  1.  The  energy  balance  at  the  surface  establishes  the  boundary  condition  for  solution 
of  the  boundary  layer  equation*.  The  energy  balance  yields. 


S  ■  -  S,  -  “vhv  *  l .  (V'v  'V1-  *  L  (V>w  (bkn)v  *  ev°Ti  *  °v  -  "rhc 

k  k  o 


vhere  the  mas*  balance  across  the  surface  la  given  by. 


(1) 


with 


a  ♦  b  •  a  ♦  m 
w  r  c  g 


(2) 


\m\.  (v’w  *  hc.  mzrl,  (v\,  (v 


•  ^  k' 


tad 


\ '  k-  ;  b«. '  r~  ("kM>v ( 


The  dependency  of  both  CH  and  CB  and  implicitly  cc  hv  Is  seen  in  the  boundary  layer 

d 

expression*.  For  the  surface  energy  balance  It  le  expedient  to  define. 


p.  ue  CH. 


(3) 
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Fl<ura  1,  flurface  Kntrgy  Balance  Mechonlene 


Cubit ltut lag  for  from  thi  boundary  layer  eapreieione,  obeervee, 


Surfed  Coii vio t Ion  eeniibie  Chemical  energy  tronefer 

conduction  energy  tr«nir«r  with  eeeciciatea  with  eurrece 

tim  CH  dependent  on  rnotloni  from  ■•••  tronefer 


lini«y  tronefer  to  Energy  IrwiTir  to  the 
thi  eurfeot  by  curtate  by  Internel  |«i 
oondeneed  material  producte 


(b) 


»•  hwD 


♦  a  6 
v4rad 


*n«r«r  trial  fir  wi y  Herodletlon  Abeorptlon  Energy  trim  fir 

froa  lurfioi  by  of  Incident.  b  T4ohenictl 

fi  ifBli  <14  IIUI  flow  rirttitlon  ii  <•  Ion 


1>n  above  equation  Include!  thi  irfiot  of  i  non-unity  lewla  number  and  li  applicable  for  both 
lialnir  end  turbulirt  flow,  Hie  energy  balanoe  ouuylea  thi  boundary  liyir  aolutlon,  which  It  required  to 
obtain  epeolee  concentratlona  and  the  tramfer  ooefflolen-e ,  to  the  lndepth  loJutloni  of  the  oharrlni 
ablator  analyili. 


(b'  :•  fineril  although,  u  preuntid,  then  iohe»l oil  equljlbrlua  li  aaeumed  In  the  boundary 
layer  and  at  the  ablatln<  lurfaoe.  However,  eurface  nan  litlon  rate,  ou  m n  be  ruction  rate  controlled 

by  reaction!  bet'"-  'he  oondenied  ablator  ohar  react  a Q.e  and  thi  boundary  layer  ipeolei.  At  1  m  aiirfaoe 
teaperaturea  r»  are  often  rate  controlled  and  (mfl  -  li  determined  froet  an  Arrhenlua  rate 

expreailon.  At  iii<lnr  temperature!  e  rate  of  Injection  of  a  Into  'he  boundary  layer  nan  be  onntrollid 

eolely  by  the  rtlffuelur.  rate  of  the  boundary  layer  reeotante  td  the  reeotln*  ablatur  aurface,  In  that 
diet,  boundary  layer  eolutlona  nuet  be  determined  with  theroochenloe)  equilibrium  reaction!  Involving  the 
•  urface  Biee  balance  apeclea,  Oaaeoua  ibletlui  product!  are  detenalned  through  the  elaulltneoue  eolullnn 
of  {•)  ill  iarfioi  th*r»oo)i#ol cil  i^ul llbi luft  :iiotlun»  ixpriMid  In  l^rvu  of  pirltAl  firtfliurtfl  ind 
equilibrium  oonetontei  (b)  equetlone  of  ooneervitlon  of  element!  epeniee  at  the  eurfaoei  end  (o)  heltone 
partial  praeeure  law. 
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Th*  equation*  for  eoneervatlon  of  elemental  tpenlte  com*  directly  fro*  th*  roniervet  1  on  of 
•pool*  equation*  end  tbo  aiiC/tlt  eurfeee  reaction  boundary  condition  (l?,  2}).  Hck'e  lew  nuet  be 
teeuaed  In  order  to  daaenetrat*  the  tut  colutlon  fora  in  the  turbulent  chi  (21).  However ,  et  ell  tine* 
the  betaideiy  layer  reaetenta  ere  aeettMd  te  fa  taidargolng  eloulteacoue  equilibrium  dlceoeletlon- 
recombination  reaction*  eaong  theneelve*  et  the  eurfeee.  In  feet,  e  continue*  between  the  two  ertremea 
of  diffueloa  Halted  and  rate  Halted  reaction*  nay  be  effeeted  ee  long  ee  the  procaaati  ere  confined  to 
the  eurfeee.  Therefore,  reeetlone  between  eoae  ablator  epeolet  end  the  beieidary  reactant*  nay  be  dlffualcn* 
Halted  while  ether*  ere  timultanaouely  rate-limited  (8,  21,  2b). 

1h*  ablation  aeehanlea*  verlouely  termed  eubliaetlon,  vepori ration  an  1/or  Halting  with 
epallatlon  ere  eharaeteriatloally  represented  ee  change  phase  reeetlone  with  energy  ebiorptlon.  Tlie 
•wbllaatlea  or  raporlietlen  le  generally  aeeoelated  with  the  eheateel  reeetlone  between  ambient  boindary 
leper  eenitltuente,  pyrolyele  geee*  end  the  eubllalng  eurfeoe  aeterlal,  Consequently,  the  heete  of  pheee 
change  ere  aoet  often  Included  with  ehcaleel  enthelplee  of  formation  of  the  eurfeee  oner  eondeneed  pheee 
epenlee,  Theee  pheee  change  reeetlone  effectively  deoreeee  the  ehealeel  reaction  enthalplee, 

Ttie  ablation  aeehanlea  termed  aeehenleel  e  roe  ion  Include*  liquid  layer  flew,  aeehanlea] 
epelletlon,  and  a*  eh  an  leal  cheer  of  eurfeee  arterial,  liquid  layer  ablation  te  dteeueeed  In  the  neat 
subeeetlca,  Mechanic*]  epelletlon  of  charring  ablator*  can  be  eeueed  by  the  internal  preeeur*  of  the 
gee  a ■  pereelatlng  through  the  char  (27),  The  etreaeee  sen  eeeaed  the  tenet le  etreagth  of  the  ohar  for 
oaay  low  char  etreagth  ablator*.  At  a  cooiequenoe,  a  part  or  aoet  of  the  char  can  eaparat*  from  th* 
virgin  aeterlal.  Meehan! rai  rliaar  ablation  la  aharaetarlted  h“  aaall  aelid  char  pertlele*  leaving  th* 
eurfeee  due  to  eurfeee  ehear  fo/ee*.  goth  of  thee*  proeeeeet  nan  b#  expedited  by  th*  Impingement  of 
eaall  eolld  pertlele*  in  th*  boundery  layer  ae  frequently  oecur*  in  th*  rocket  nettle  eroelon.  Mathanloal 
cheer  abletien  1*  normally  aeeoelated  with  mart  eever*  enrlronaental  ronditlone  end  thinner  eher  layer*. 

Kq  (b)  1*  not  readily  rtoognliabi*  In  the  fora  n*  preeented  but  ean  b*  rearranged  to  an 
alternate  fora.  An  alternate  oholc*  for  th*  eurfeee  energy  balance  1*  th*  following  expreaelen  (19,  2f>)i 


S  •  p.u.ldf)  -  <h,)  1  -  *  Vrt4  *  ioh,al**l  *  A,«  *  V»0  *  ^blocking  (J) 


fur*  convection  tern  Mae*  tranefer,  vaporleatton  bicoking 

with  Cu  bearing  no  end  eurfeoe  eheaietry  tern  term 

dependence  on  *n  or  h^ 

where 


Th*  b'eehing  term  trie**  froa  eepreeeing  i  lu  term*  of  C..  lneteed  of  0  .  Where  c((  1*  th* 

.  0  0 

Blanton  number,  0H,  with  taro  aaaa  Injeetion  (i.a,,  «v  ■  0).  Thl*  rellevet  the  reepteielbi  llty  of 

determlalag  th*  pur*  eaaveetlv*  Itanton  i  ebar  from  eeJutiona  to  tpeelae  coupled  eonewrvetion  aquation*, 
tapir 1**1  re  1  at  lane  frequently  feel  111  a* «  determination  of  th*  blocking  term. 


in  neat  *****  tv,  a,  and  ar*  eonetanta  while  q^U)  le  a  table  arriving  from  th*  chock 

layer  tolutlent.  At  elated  earlier,  quit*  frequently  q  li  evaluated  by  Eckert'*  reference  enthalpy 

0 

nethod  for  t  few  elapl*  ehepce,  both  for  turbulent  and  laalnar  enrlronannta  (19,  t(>) . 
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2.3  Ablation  with  Liquid  Layar  formation 

Thl*  Mehulw  may  b«  considered  at  «n  Intermediate  prootit  to  tnt  ablation  vlth  only  turfaca 
mass  tract  far  and  tha  aora  ooaplax  process  Involving  daeoaipotltion  In  depth.  Tha  eaaa  of  ample  salting 
with  tba  liquid  layer  initantanaoutly  avapt  fros  tha  surface  or  total  vaporisation  at  tha  llquid-ga* 
Interface  for  vary  high  vlaooalty  satarlala  say  ba  treated  by  tha  sethoda  of  tha  previous  aaotlon.  How¬ 
ever,  oaaaa  of  ooabinad  salting  and  evaporation  of  glaaay  satarlala  vhara  tha  liquid  layer  partially 
float  around  tba  body  and  partially  evaporates  require  a  sort  ooaplax  sodel  of  thle  naohanlaa.  Exact 
solutions  for  tha  tranalant  atata  arbitrary  heat  flux  ara  given  In  (20,  29,  JO)  at  stagnation  aa  vail  at 
around  tha  body  for  arbitrary  valoolty  and  taaparatura  gradient  In  tne  liquid  layar.  Thlo  aolutlon  la 
eoupled  vlth  Internal  radiation  tranaport  In  (31)  and  a  revlev  of  the  problem  la  given  In  (12). 

for  epeolfle  retulto  of  tha  liquid  layar  anaiyala  and  bibliography  tha  reader  la  directed  to 
(20,  29,  30,  31 1  3b).  It  should  be  noted  that  tha  aathoda  reoosnended  ara  applicable  to  any  material* 
displaying  liquid  layar  flov  around  tha  body  and  not  Juet  gli iay  aaturlala.  Tha  isportanoa  of  tha 
Internal  taaparatura  dlatrlbution  la  partleularly  noted  baeauaa  of  tha  atrong  affect  of  taaparatura  on 
vlaooalty  of  tha  liquid  layar  and  thue  on  ablation  ratoa. 


It  la  in  thle  area  that  technology  haa  exerted  ooneidereble  effort  to  arrive  at  phyaloally 
meaningful  and  raaeonably  acourate  aolutlona  (0,  19,  21,  2b,  33). 

Actually,  tha  aurfaoa  aquatlona  and  tha  oouplad  boundary  layar  aolutlona  naraly  rapraaant  a 
eat  of  boundary  eonditione  aattlng  tha  ataga  for  tha  ooaplax  tranalant  parabolic  heat  conduction  aquation 
daaorlblng  tha  indapth  phanosana,  that  la, 


Tha  usual  aaausptiona  constraining  tna  Indepth  solution  ara  aa  follovsi 

1.  On#  continuous  aquation  applies  aeroaa  tha  entire  ablator  to  lnoluda  tha  deposition, 
porous  ohar,  decomposition  (or  reaction)  and  virgin  materiel  ragloua. 

2.  Tlia  pyrolyflng  material  either  daooapoaaa  oosplataly  Into  gaseous  apaclaa  or  laavaa 
partially  a  solid  residua  capable  of  malting  or  entering  Into  heterogenous  rsaotlona  st 
tha  surface. 

3.  Tha  Mae  flux  of  indaptb  pyrolysis  gaseous  products  ara  aanaad  Instantaneously  at  tha 
aurfass  after  deposition  and/or  decomposition  reactions. 

b.  The  pyrolysis  gaaaa  ara  alvaya  In  thermal  aquillbrius  with  the  porous  char. 

3.  The  density  of  tha  pyrolysis  gases  are  aseuaeo  email  In  comparison  to  the  density  of  the 
condensed  solid  phase. 

6.  Pyrolysis  gas  enthalpies  a re  considered  known  functions  or  looal  temperature  and  external 
pressure. 

7,  The  thermal  conductivity  and  epeolfle  heat  of  the  eoitde  can  be  represented  as  linear 
funotion*  of  the  condensed  phis*  denelty  and  functions  of  taaparatura. 

0.  The  ablation  acdel  aasuaee  tha  one  dimensional  tranalant  heat  and  mass  tlov  to  be  normal 
*o  tha  eurfaoe. 

The  velocity  dependent  term*  such  ••  tnose  due  to  viscous  stress,  friction,  mixing,  ate., 
are  negleoted  In  the  transient  Mat  oond"."tlon  aolutlon. 

10.  The  ablating  surface  will  serve  as  the  reference  frame  for  the  moving  coordinate  system. 

11.  Thermal  deeoeipoeitlon  and  deposition  are  assumed  to  follow  eimpl*  Arrhenius'  rata 
expressions. 

12.  Internal  radiation  hast  absorption  in  a  aeml-trareparent  allleloue  char  la  usually 
represented  ae  a  gradient  expreiilon  added  to  the  heat  transfer  equation. 

The  preceding  are  generally  ecoepted  aa  tha  standard  assumptions  accompanying  an  indepth 
ablation  aolutlon.  hovevar,  other  options  exist  presenting  Internal  pressure  distribution  effects;  trans¬ 
piration  cooling-,  and  model*  for  evaluating  tha  internal  radiation  contributions. 

Upon  exercising  the  above  assumptions  one  may  obtain  directly  the  following  transient  Indepth 
reacting  host  conduction  equation  In  one  dimensional  partial  differentia)  equation  form  (;’l,  3b)  i 
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with  the  auxiliary  equations 

3p  3m 

TTmJ?  *«  °dp*0c’°c  ll0) 

Enthalpy  of  decomposition,  .  and  enthalpy  of  ueposltlon,  Ah^ ,  are  f(P^,T).  Typically,  an 
example  decomposition  reaction  might  be  1 3,  jb). 


r<ta  Ni  Fc  1  *  U'D(E%  ?c  )  AHdfi  sfTlR^fCDEP^glTlE^^KDP 
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Virgin  Composite 


Pyrolysis  gaseous  products 


eould  char 

such  that  Pdf>  «  *  £>li;)  *  (l-Do2  and  T  being  the  volume  fraction  where  c>i2  and  o?  are  completely 

converted  into  gaseous  products. 

Then  the  governing  Arrhenius  rate  expression  might  be,  for  Instance, 
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•  *"  *  u-r»k2n2(0) 


where  f  ■  - 
o  o. 


virgin 


nimllarly  ^or  a  typical  deposition  reaction  which  might  be  2E,,R2  (pyrolyele  gae)  °^c  EP^ 
(further  pyrolyele  gas)  ♦  3E  (solid  char)  the  governing  Arrhenius  rate  expression  would  probably  be 


— «  -k  ft  (0) 
ot  c  c 


-  p'-  n  p  (t)  "c 

>  •  f  (pj  c  !^rl 


Note  the  arrow  in  the  above  deposition  reaction  shows  going  to  completion  as  stated  In  the 
approxlmetloni .  A  reverse  reection  would  involve  equilibrium  conetents  and  partial  pressure  dependence. 

Considering  this  es  the  complete  picture  then  conservation  requires,  as  seen  by  the  auxiliary 
equations ,  that 


cm  3p  3d.,  3p 

_&  .  _•  .  -M  ♦  — £. 

3y  3t  3t.  3t 


(14) 
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finally  the  char  transparency  radiation  tar*  aa  proposed  In  (2)  would  be  quite  adequate  for 
this  typical  analysis,  namely 


3F 

3y 


a(J-J  )  ♦  e( l-p  )F 
e  o  a 


#-a(y-<j)  +  R  #-a(21-y-c) 


o 

»2A(X-o) 


where  J  Is  solved  from  the  auxiliary  set  of  equations 


afi  .  2(j-j.) 

2  *  e 

3yz 


2J 

3y 


(15) 


(16) 


(17) 


J 

e 


Uu^oT 


1> 


with  boundary  conditions  J  (y»o)  »  a^q^,^ 
and 

F1 

and 


«^V  • 

"  a  1+R  “wS-ad 
o 


(18) 

(19) 

(20) 


J(y-X)  -  fjU)  -  0  (21) 

All  of  these  equations  are  thence  assembled  in  finite  difference  fora  and  solved  numerically 
on  the  computer  (21,  36).  The  solutions  provide  accurate  theraal  histories  and  profiles  throughout  the 
ablator  as  veil  aa  a  realistic  prediction  of  the  decomposition  and  deposition  process. 

2.5  Mechanical  Removal 

It  Is  known  that  some  ablation  materials  under  certain  external  or  internal  hlgh-preaau  i 
conditions  may  either  experience  mechanical  opallatlon  of  the  surface  or  other  failures  resulting  In 
mechanical  removal  of  nterlai.  The  mechanisms  or  failure  criteria  for  such  ablation  processes  have  been 
studied  (27,  37),  criteria  assumed  and  postulated,  and  their  results  Incorporated  In  the  thermal  ablation 
models,  further  detail  will  be  provided  In  Section  3.0. 

2 . 6  Specialised  Ablation  Models 

The  sample  model  of  ablation  material  exhibited  earlier  to  represent  the  typical  lndeptb 
decomposition  end  -  -position  mechanisms  vas  symbolically  that  or  phsnollc  nylon  In  absencs  of  mschanlcal 
removal,  The  symbolic  pyrolysis  gases  rsleased  by  decomposition  of  this  model  were  a  four  component 
mixture,  namely  KR^,  E^g,  Rg,  and  P  representing  respectively  methane  (CH^),  acetylene  (C^),  hydrogen 

(Hj) ,  and  a  combination  of  assumed  nonparticipating  gasea  (l.e. ,  non-reacting  with  boundary  layer  species) 

such  aa  nltrogsn  and  hydrogen  cyanide  (35).  In  this  model  the  completely  volatile  nylon  reinforcement  vss 
represented  by  0g  whereas  the  (completely  volatile)  and  combination  represented  the  density  of  the 

phenolic  resin  with  T  being  the  volume  fraction  of  phenolic  resin  In  the  entire  virgin  plastic  composite 
(3,  21,  3b).  The  last  reference  treats  a  rather  sophisticated  coupled  chemically  reacting  boundary  layer/ 
charring  ablator  problem  vlth  Interesting  results,  but  unfortunately  not  quite  tractable  by  a  practicing 
engineer.  With  the  equilibrium  surface  reactions  between  these  pyrolysis  gaseous  products  and  the 
boundary  layer  apecles  In  mind,  tbe  choice  of  vhlch  components  era  representative  or  most  significant 
becomes  largely  an  art  depending  on  the  range  of  surface  temperatures  (It).  The  problem  of  handling 
coneecut lve  equilibrium  surface  reaction*  Involve*  no  more  than  writing  one  combined  stoichiometric  ex¬ 
pression  and  the  associated  effective  equilibrium  constant.  0c  the  other  hand,  simultaneous  equilibrium 
surface  reactions  momentarily  appearing  aa  primary,  secondary  and  tertiary  In  some  order,  only  to  have 
their  roles  reverted  or  permuted  at  the  temperature  changes,  presents  the  only  real  challenge  for  choice 
of  original  pyrolysis  components.  This  certainly  vould  seem  to  preclude  the  four  component  model  from 
even  representing  the  ablation  of  phenolic  nylon  much  lets  the  general  ablator. 

However ,  knowledge  of  the  thermoc hem leal  properties  of  materials  such  e*  teflon  or  carbon  and 
the  graphites  and  tbtlr  respective  ablation  behavior,  in  such  veritable  detail  permits  specialised 
handling.  This  greatly  simplifies  the  ablstlon  equations  lurther  facilitating  Inclusion  of  more  detailed 
chemistry  (39).  For  instance,  carbon  (graphite)  can  form  the  following  ablation  products  limited  in  the 
okldatlon  and  sublimation  regimes:  (b,  bO,  bl,  US)  c,  Cg,  CJ(  C^,  C^,  Cp,  CQ  ,  C^,  and  CjJIj.  Kor  qusrtl 

the  chemistry  Is  much  simpler  as  the  only  reaction  requiring  consideration  1*  (2b). 


lit 


Si02(t)  l  SlO^Og 

where  (1)  Aerate*  liquid. 

Intermediate  etepe  can  include: 

Si02(t)  l  SiOg 

SiOg  *  SKOg 

sio2  l  SiO*|o? 

For  graphite  vhere  the  eurfe.ce  reaction*  and  hence  ablation  ia  reaction  rate  Halted  the 
following  heterogenous  completion  reactions  are  representative,  namely: 

C  (a)  ♦  0  ■»  CO 

Simultaneous  completion 
(reaction  rate  limited) 

C  (s)  ♦  |  °2  *  CO 

The  mass  fl’tt  of  C  (*)  leaving  the  surface  by  interaction  with  0  and  0?  being  represented  by 
the  combi:  d  rate  law  and  Arrhenius  expression: 


ni  «c  "2  flc  iff 

a  •  const  P  1  (—  (C  )  ]  [—  (C  )  )  exp*1  (22) 

fl0  M  Ro2  2W 

Similarly  for  teflon  vhose  ablation  and  surface  erosion  occurs  at  the  surface,  the  folloving 
rate  dependent  formation  of  the  teflon  monomer  and  a  subsequent  oxidation  reaction  can  be  expressed 
a*  (U3). 

(C  F  )  (s)  -  N(CF  )  (s)  l  N(CF  )  (?)  'I  2N(CF  )(g) 

N  2  ‘  2 

K,„  (23) 


(CFg)(g)  *  0  *  COFg 


(vith  arrow*  in  both  directions  for  diffusion  limited  and  Harvard  direction  only  for  rate  limited 
reaction!) . 


The  basic  advantage  of  considering  the  teflon,  o.uartz,  and  carbon  models  separately  from  the 
generalized  model  lies  in  the  feature  of  having  to  consider  only  Indepth  transient  heat  conduction 
Internally.  All  reactions  are  confined  to  the  surface  area. 

The  silicone  elastomers  conversely  pursue  the  opposite  route.  Ablation  occurs  primarily 
through  the  indepth  decomposition  mechanisms  vith  relatively  sparse  surface  eroelon  reactions.  The  porous 
silica  char  produced  provides  a  highly  lneulative  layer  that  primarily  reradiatec  energy  back  to  the 
environmental  flov  field  (Ul»),  The  basic  indepth  reactions  consist  of  unzipping  of  the  highly  cross- 
linked  bonds  of  the  elastomer  vith  the  release  of  a  fev  basic  pyrolysis  gases  such  as  water  vapor,  carbon 
monoxide  and  methane  (T).  A  highly  pc roue  silica  char,  SiOg  nay  reduce  to  silica  carbide  and  carbon 

monoxide  later  capable  of  entering  into  reaction  vith  the  boundary  layer  constituents  end  forming  the 
dioxide  again  (7.  1*5,  b6). 

The  silica-resin  composites  present  the  combined  complexities  of  the  silicone  elastomeric 
indepth  reactions  with  the  phenolic-nylon's  phenolic  resir.  decomposition  ( U6 )  and  quartz's  surface 
reaction  characteristics  (2b). 


The  carbon-graphite  resin  composites  essentially  proceed  vith  the  same  or  comparable  indepth 
reaction  kinetics  as  the  phenolic  nylon  model  vith  char  erosion  surface  reactions  similar  to  the  pure 
graphite  model  (39,  b7). 


Radiating  Shock  Layers 


At  entry  velocities  above  b^.OOO  fps ,  stagnation  point  heating  can  be  primarily  caused  by 
strong  atomic  line  and  continuum  radiation.  Under  these  conditions,  ablation  rates  result  in  msss  fluxes 
from  the  heat  shield  on  the  order  of  20  percent  of  the  air  flux  through  the  shock  (b8,  kg).  Classically, 
such  conditions  result  in  "blovn-off"  shock  layers  with  the  convective  heating  reduced  many  orders  of 
magnitude  below  the  value  that  vou;d  occur  without  bloving. 


Since  most  of  the  coupling  radiation  occurs  in  the  ultraviolet  spectrum,  ablation-products  in¬ 
jection  will  not  offer  any  significant  alteration  of  these  phenom  na,  although  they  will  be  somewhat 
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accentuated.  These  conclusions  we'e  borne  out  by  Coleman ,  Lefferdo,  Hearne,  and  Vojvodleh  (18)  where 
radiative  coupling  to  ablating  heat  shields  were  studied ;  reductions  In  radiative  heating  from  ablation- 
products  infection  for  the  cases  considered  by  Rlgdon,  Dlrling,  and  Thomas  (1*9)  reach  30  percent. 

The  recent  decrease  In  emphasis  on  msaned  entry  missions  at  velocities  in  excess  of  36,000  fpe 
have  witnessed  a  corresponding  decrease  in  the  attempts  to  further  understand  and  analytically  model  the 
phenomena. 

3.0  SYNERGISTIC  ABLATION  MECHANISMS 

It  Is  known  that  seme  ablation  materials  under  eertaln  external  or  internal  high  pressure  con¬ 
ditions  may  either  experience  mechanical  spallation  of  the  surface  material  strata  or  other  failures 
resulting  in  mechanical  removal  of  material.  The  mechanisms  or  failure  criteria  for  such  ablation  pro¬ 
cesses  have  boen  studied,  (i>7,  37)  criteria  assumed  and  postulated,  and  their  results  Incorporated  in  the 
thermal  ablation  models.  Although  these  materials  have  been  studied  analytically  with  some  vigor,  recently 
developed  blgn  pressure  test  facilities  have  revealed  a  large  deficiency  In  analytical  prediction*  of  the 
obse.-ved  mass-loss  data.  As  a  consequence,  analysts  have  resorted  to  empirical  and  semi-eoplrlcal 
correlations  of  the  test  data  with  rather  questionable  results  when  extrapolated  to  actual  flight  environ¬ 
ments  . 


Understandably,  empirical  correlations  of  test  data  are  at  best  provisional;  little  Insight 
into  the  actual  mechanisms  operative  during  the  ablation  processes  ia  gained  (50,  51).  Conversely,  rather 
sophisticated  analyses  of  the  thermochemical  erosion  of  a  variety  of  carbon-based  materials  has  evolved 
(3,  52,  53).  Development  of  these  analyses  las  been  closely  followed  by  a  number  of  papers  dealing  with 
various  aspects  of  the  thermomechanical  ablation  problem  (37,  27.  35)  and  are  specialized  to  the  oarticular 
material  unaer  study.  The  composite  of  these  analyses  has  failed  to  provide  repetitive  correlation  of  the 
test  data  —  obviating  their  generality.  It  Is  now  evident  that  a  large  number  of  disciplines  are  required 
to  effectively  approach  the  analytical  description  of  even  what  previously  was  described  as  homogeneous 
materials,  graohlte. 

Graphite  response  mechanisms  will  be  used  here  as  a  typical  example.  Analogous  response  would 
occ-r  in  some  other  materials. 

3 . 1  Environmental-Ablation  Interaction 

3.1.1  Typical  Nose  Tip  Interaction 

The  definition  of  material  response  that  1#  required  to  account  appropriately  for  the  Inter¬ 
action  of  the  local  aerodynamic  environment  at  the  exposed  material  surface  and  for  the  resulting  lndepth 
phenomena  1*  complex.  To  illustrate  this  behavior  In  an  applied  situation,  a  r.sce  tic  .3  analyzed  for  a 
ballistic  entry  flight  case  perturbing  first  the  extent,  of  mechanical  erosion  and  then  transition 
criteria.  The  results  are  presented  In  Figure  2. 

The  results  of  the  sensitivity  analysis  for  the  cases  examined  are: 

Case  1.  The  baala  for  the  noae  shape  change  results  in  this  flight  simulation  case  la  a  sil.  coua 
note  tip  material.  A  nominal  transition  criterion  neglecting  roughness  and  blowing 
was  employed  to  determine  transition  to  turbulent  heating.  Mechanical  erosion  waa 
deduced  from  ground  test  data  report  in  (51*). 


Case  2.  This  case  shows  the  sensitivity  of  mechanical  erosion  effects  on  nose  shape  change  (55) 
and  recession.  The  case  Is  Identical  to  Case  1  with  the  exception  that  the  mechanical 
correlation  utilized  vas  for  a  factor  of  3  increase  grester  In  recession  over  that 
employed  with  the  sllaceout  material  In  Case  1  at  a  stagnation  pressure  of  100  atm. 

Case  3,  This  case  retains  all  assumptions  of  mechanical  erosion  for  the  silaceous  material 
used  in  Cass  2.  The  only  exception  Is  that  laminar  flov  la  assumed  to  persist 
throughout  the  reentry  flight. 

Case  U.  This  case  employs  the  same  material  assumptions  as  Cases  2  and  3,  with  the  exception  that 
turbulent  flow  is  presumed  to  exist  from  the  start  of  reentry  throughout  the  flight  period. 
The  all-turbulent  run  indicated  here  simulated  a  fully  turbulent  stagnation  point  through¬ 
out  the  flight. 

In  summary,  the  results  reflected  in  Figure  2  show  large  changes  in  spatial  and  temporal  nose 
shape  configurations  due  to  the  tbermomechanical  erosion  assumptions.  Extremes  In  all-laminar  and  all- 
turbulent  boundary  flov  produce  signlflcan.  effects  upon  both  the  final  shape  and  total  tip  recession. 

When  using  the  some  r.echanlcal  erosion  correlation,  the  nominal  transition  criteria  naturally  lie  between 
these  two  extrsmes  and  are  more  closely  represented  by  an  all-turbulent  layer.  The  significance  of 
mechanical  erosion  and  flow  field  assumptions  are  clearly  evidenced  by  this  Illustrative  case.  For  purposes 
of  this  paper,  primary  attention  la  focused  on  mechsnlcsl  erosion  and  the  relatively  nev  approaches  used 
to  analyze  and  gain  Insight  Into  the  mechanisms  promoting  this  erosion. 

3.1.2  Afterbody  Response 

Recently  differential  ablation  patterns  have  been  observed  or.  afterbody  surfaces.  Previously, 
intriguing  ablstion  patterns  on  meteorites  and  tsktltes  had  bean  noted  and  studied  (e.g..  References  56, 

57,  58).  Later,  attention  was  focused  on  cross-hatching  and  related  roll  torques  by  the  advent  of  small 
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hyparveloc lty  vehlclai  uj  their  erratic  roll  behavior.  Larson  and  Hateer  noted  the  development  of  croea- 
hatohlng  oo  camphor  aodels  Into  a  pattern  resambllng  that  obaerved  on  many  aeteorltea  ($9).  Canning 
Wilkins  and  Tauber  have  reported  on  turbulent  vedgee,  streams lie  grooves,  and  cross-hatched  pattern*  fro* 
their  balllatlc  range  and  vlnd-tunnel  teata  and  propoaed  a  phyaioal  flov  model  to  explain  the**  pattern* 
(60.  61.  68). 

Laganelli  and  Haatlar  have  reported  the  results  of  teete  In  the  Malta  rocket  exhaust  facility 
and  tha  Langley  8 -ft  structure*  tunntl  on  a  vlda  varlaty  of  aatarlals.  ranging  fro*  wood  to  carbon 
phenolic  (63).  Particular  aaphasls  vaa  placed  on  char-forming  materials. 

williams  (68)  (Figure  3)  obtained  croaa-hatchlng  In  a  hypersonic  vlnd-tunnal  on  aaveral  low 
tMperature  ablators,  mostly  subliming  camphor  and  Korotharm.  H*  shoved  th*  seal*  of  croaa-hatchlng 
dacreaalng  with  Increasing  prsaaure  for  camphor  and  the  pattern  of  melting- vaporising  naphthalene  to  be 
similar  to  that  of  subliming  camphor,  A  uni qua  feature  of  the**  tests  vas  the  direct  measurement  of  roll 
torques  during  the  ablation  procssa.  Later  teata  (6$)  extended  this  work  to  a  more  realistic  simulation 
of  the  reentry  vehicle  flight  regime  and  geometry.  Similar  results  vers  obtalnvd  in  camphor  and  vax  by 
Stock  and  Olnoux  (66).  NcDevitt  (67)  Investigated  the  coupling  betveen  roll  dynamics  and  ablation  pattsrns 
by  testing  spinning  Korotharm  and  ammonium  chloride  models  in  e  hypersonic  vlnd  tunnel. 


Figure  3.  Croaa-hatchlng  on  a  80°  Camphor  Cone  at  Mach  6 
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Various  analyses  hat*  basn  uli  to  explain  tho  phsnertenoo  of  cross-hatching  i  only  t  f*v  cob  bo 
mentioned  boro.  Tobol  (68)  auggeated  vortteaa  u  tho  origin  of  ei'oaa-hetehlc*.  lechtabelm  (69)  bellavae 
o  thin  liquid  layer  It  oooootlol.  Probsteln,  Oold  end  Soulleo  (70)  believe  that  crust— hot  shine  oceuro 
boo oust  of  *  resonance  between  boundary  layer  disturbance*  and  a  suit obi#  Inelastic  dtforeabl*  surface . 
la  ob  in  too  t  loot  loo  of  tho  pressure  and  temperature  distributions  over  a  tvo-direns  local  non  ablating  vary 
wall.  Willi aaa  and  Infer  (71)  bar#  theoretically  proved  and  experimentally  confirmed  a  vail  toaporature 
phaala*  eorreiboc donee  to  boundary  layer  preaaure  disturbances.  Williams  and  Infer  (71)  also  report  that 
a  pre-eroaa-hatched  nocablatlnf  Model  conflnsed  cross-hatching  experience  by  tndlcatlnf  higher  heat  trena- 
fbr  to  the  vallqre  than  to  the  ridges t  this  la  diametrically  opposed  to  normal  supersonic  heat  transfer 
experience  which  shove  peek  heating  an  protuberances. 

3.2  Oranhlts  Eros loo  Mechanisms 

It  has  been  observed  by  a  number  of  investigators  (72 ,  73,  7b),  that  \ader  certain  prassure- 
teidperature  anvironaanta  bulk  polycryatalllne  graphite  undergoes  loss  of  solid  aetorlel  of  minute  quantity 
characteristic  of  the  constituent  else  of  tha  notarial  components.  In  o  closer  examination  of  tho  type  of 
behavior  that  occurs  vhan  graphites  are  subjected  to  high- temperature  environments,  a  cross  section  of  a 
tested  specimen  Is  shown  in  Figure  b.  The  figure  reflects  an  lndepth  degradation  of  aetorlel,  shoving  a 
rather  open  Irregular  structure  proceeding  into  tha  eat trial  mud  the  appearance  of  an  interface  shoving  a 
lata  porous  structure  with  Increased  material  aasa.  A  rlev  of  tha  surface  demonstrates  large  difference# 
In  materiel  structure  and  preferential  lota  of  material  revealing  irregular  rectangular  appearing 
particles. 
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Figure  b. 


high  Temperature  Exposed  Cross  Section  and  Surface  Photomicrograph  of  Graphite  (f'TJ-S) 


Scxae  perspective  can  be  focused  on  the  results  of  the  ablated  pattern  illustrated  in 
figure  W.  Graphite,  by  its  vary  nature.  Is  an  inhomogeneous  oa* trial  composite  consisting  of  filler 
particle#  and  binder  structure  vith  a  high  degree  of  porosity.  The  differential  nature  of  the  constituent 
Materials  used  in  fabrication  of  tha  graphite  Is  believed  to  result  in  preferential  erosion  giving  rise  to 
an  irregular,  roughened  surface  that  Interacts  closely  to  the  environment  at  the  heated  surface  of  the 
material.  This  interaction  involves  a  couple*  association  of  the  local  environment  and  the  material 
behavior. 


3.2,1  Physical  Characteristics 

An  understanding  of  tha  microstructure  of  graphite  requires  sees  knowledge  of  the  processes  by 
vhlch  graphite  it  formed.  The  conmerclel  graphite*  used  in  aerospace  applications  can  be  described  as 
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synthetic  or  nsnufacturad.  These  bulk  graphites  oro  manufactured  by  mixing  am  ergenlo  binder ,  usually 
o  cool  tor  pit  eh  |  with  earoon  filler  particles  and  aoldlag  or  extruding  the  fixture  to  fora  o  "green* 
billet.  The  hinder  la  then  carbooised  bgr  hooting  to  te^erotureo  on  the  order  of  1300T,  followed  by 
further  hooting  to  temperatures  up  to  JkOO’F,  to  preduee  o  no  re  ordered  aleroetruoture.  The  filler 
aoterlol  nsy  he  any  portieulote  graphite  or  eorhon,  eueh  oe  ooturol  graphite,  carton  blook,  coke,  etc. 
for  eeenonie  reoeone,  filler  aoterlol  la  noet  bum only  node  fron  the  coked  realdue  of  the  petroleum 
refining  prooeet.  The  cote  note  le  crushed  tad  grow*  to  eliniaote  greet  poroolty,  end  tereenod  into 
controlled  particle  elxei.  Because  of  unavoidable  lntercryetalllne  mieroports,  the  particles  retain  . 
minimum  poroolty  of  5  to  8  percent  (75),  reeulting  in  apparent  demit  lee  of  about  2.0  to  2.2  gn/cc. 

Figure  5  provide*  o  echenatle  of  the  propoeod  hulk  graphite  arrangement  after  grtphltleotloe 
(the  dneerlptioo  given  here  abet*  the  excellent  dleeueeloa  of  graphite  by  Filey  (78)).  The  filler 
particle*  appear  os  an  array  of  loosely  pocked  irregular-shaped  coke  elenent*  can  anted  together  by  the 
carbonised  binder.  The  ratio  of  hlnder-to-flller  weight  nay  be  varied  over  wide  Unite,  with  nornal 
practice  being  between  20  to  l»0  percent.  In  the  baking  cycle  (l$00eF)  the  pitch  binder  is  pyrolyscd, 
losing  about  bO  to  SO  percent  of  its  initial. weight.  During  tbs  baking  cycle,  end  upexi  subsequent  grapbl- 
tlsotlon,  the  hinder  shrinks.  Since  shrinkage  it  constrained  by  the  filler  particle*,  the  volunt  nonin ally 
occupied  by  the  hinder  hea  a  large  void  content.  The  resulting  nacroporea  in  tbs  binder  allow  for  oaay 
passage  of  evolved  gases  through  the  hinder  natrlx  and  account  for  the  reletlvoly  high  porosity  of  hulk 
graphite.  Each  filler  particle  la  held  in  the  reeulting  natrlx  by  a  network  or  binder  bridge*  occasion  ally 
interrupted  by  interconnected  necroporet.  The  nacropo re  dleneter  varies,  but  is  about  the  sans  as  the 
disaster  of  the  filler  particles.  A  lias  of  Interconnected  sac  roper* s  nay  propagate  a  aechaalcal  failure 
since  this  le  the  path  of  nininua  cross  section.  Tbs  nicrostructur*  evidenced  in  Figure  S  erpoaei  the 
beterogeneoua  qualities  of  graphite  In  detail.  Filler  particles  consist  of  closely  packed  crystallltos 
and  lntercryntalline  nlcropores.  Significantly,  the  crystallites  are  cot  completely  ordered i  the  binder 
crystallites  arc  even  less  ordered,  vhich  account*  tor  the  nininal  anisotropy  exhibited  In  the  thermo- 
phyelcal  properties  of  hulk  graphite. 


Figure  5.  Schematic  of  Graphite  Microstructure 


3.2.2  Praferentlal^roslon 

Father  extensive  Investigations  of  the  preferential  er  sicn  aspecte  of  graphite  have  occurred 
In  recent  literature.  Kasha  and  Gchryer  (72)  observed  loee  of  particle  solid  of  graphite  nose  tip 


utiriiU  in  ««lr«nnti  typified  by  model  stagnation  pressures  of  k  to  5  eta  (taxation  pressure,  Nora 
recently ,  Made  11  (73)  photographically  thoatd  proof  of  loaa  of  particulate  aolld  graphite  in  model  itif- 
nat toe  pressures  u  low  aa  1  at*.  These  observations  hay*  not  bean  eoafinad  to  air  arc  tasting  aloes ,  aa 
evidenced  toy  the  observation  of  Wblttakar  (7*>)  vho  observed  particulate  aolld  loaa  and  e«p*rlmntelly 
trapped  the  oaolrad  aolld  aaai . 

Recant  experimental  and  analytical  laterpratatlona  of  the  prafarantial  aroaico  of  graphitic 
■atari ala  bar*  bean  offered  by  Kratach  (77)  and  McVey  {78).  Both  addraaaed  the  high  praaaura  perfomanea 
of  graphite  vhara  large  daviatlcna  in  tberaocheaical  ablation  calculation  ware  noted.  Typical  experimental 
evidence  of  the  departure  from  thermoehealeal  eroaloo  given  in  (77)  ie  ihown  in  figure  6.  Denmen  (79) 
raviewad  tit  graphite  ablation  problem  from  an  environmental  approach,  calculating  ahgpa  change  and  reces- 
al on  rataa  for  anooth  and  roughaned  valla.  It  aufflcaa  to  atate  that  the  eoopoalta  of  theae  analyte a 
remain  eemiemplrlcal  and  provide  little  insight  into  the  relation  of  the  obaerved  a roe  ion  to  the  material 
composition . 


Behavior  of  the  graphite  ahovc  in  figure  k  demon at rataa  a  complex  Interaction  of  materlalr 
composition  and  local  aerodynamic  anvironnent.  To  identify  the  materia*  quantities  that  have  an  important 
role  in  the  ablation  performance  requires  a  systematic  approach  to  laolatlon  of  individual  a  rot  ion 
mechanisms.  Tbs  approach  adopted  to  elucidate  a roe  ion  sechanlama  for  graphite  matarlala  le  schematically 
repreaeeted  in  figure  7.  The  figure  represent*  the  compositional  structure  of  graphite  vhao  subjected  to 
hyperthermal  reentry  environments. 


figure  7 


Klcrcmecbanical  Bros  loo 


In  air  tnvlrocmants,  the  boundary  layer  constituents  chemically  react  with  the  carbon  structure  resulting 
In  thermochetUcal  e  roe  Ion  by  oxidation  and  eubliBatioe  of  the  carbon  atructure,  The  poroaity  inherent  In 
graphite  afforda  a  naane  of  eubaurfaca  chemical  eroalon  to  oocur.  These  effects  encompass  subsurface 
oxidation  and  aubllnatlon  that  are  closely  related  to  the  material's  macro  and  mlcroporcsity.  Porosity 
alone  la  not  oonflned  to  that  open  to  attack  by  the  boundary  layer  constituents  but  Includes  closed  pores 
which  may ■  at  sufficiently  elevated  temperatures,  act  as  miniature  pressure  vessels  as  sublimation  occurs 
within  the  enclosed  Internal  pore  surface  area. 

Conaaquently,  the  differential  microstructure  between  filler  particles  and  binder  material 
holding  the  particles,  through  a  labyrinth  of  binder  walla  connected  and  interlaced  throughout  the  pores 
may  have  significant  effect  on  the  microoechanlcal  erosion  observed  by  graphitic  materials.  The  prefer¬ 
ential  removal  of  either  binder  or  filler  particle  material  leads  to  a  differential  recession  rate  causing 
an  irregular  surface,  promoting  eurface  roughness  and  increased  heat  transfer.  Material  protruding  about 
an  aerodynamic  atirfach  la  subject  to  aerodynamic  shear  and  pressure  forces  which,  depending  on  the  material 
strength,  extent  of  exposure,  and  bonding,  may  be  sufficient  to  erode  solid  from  the  surface.  If  the 
material  pore  structure  is  (mail  and  constricted,  a  pressure  differential  may  exist  across  the  material 
causing  local  failure  of  the  partially  eroded  material. 

To  quantise  material  composition  end  variables  which  may  be  significant  in  isolating  individual 
erosion  mechanisms  contributing  to  the  total  performance  of  graphite,  necessitates  on  approach  which  first 
dissociates  siaterlal  variables  interactions  from  the  applied  environment.  With  the  material  quantities 
characterised,  euch  as  binder  content,  filler  particle  size  and  distribution,  pore  distribution  and  shape, 
etc.,  concise  statement*  can  be  made  about  the  material  variables  influences  in  aerodynamic  environments. 
These  material  variables  can  then  be  related  to  the  observed  erosion  occurring  snd  developmental  goals 
for  improved  graphitic  materials  stated. 


3.2. j  Surface  Roughness  Effect  on  Heat  Transfer 

The  material  characteristics  are  though  to  directly  influence  heat  transfer  to  the  material 
ablating  surface.  Surface  roughness  elements  can  ocr rr  from  material  impurities  and  preferential  erosion; 
or  they  can  be  inherent  in  the  material  finish,  even  when  highly  polished.  The  effects  of  surface  rough¬ 
ness  on  heat  transfer  have  been  treated  by  numerous  investigators  (80  through  93)  with  varying  objectives. 
The  bulk  of  the  reported  studies  have  been  concerned  either  with  singular  roughness  elements  (both  two- 
end  three-dimensional)  on  symmetrical  bodies  or  with  roughness  elements  on  flat  plates.  The  utility  of 
these  data  in  a  high-pressure  environment,  and  as  related  to  heat  transfer,  has  yet  to  be  determined. 

The  data  of  Bertram,  et  al.  (91,  92,  and  93)  would  appear  to  be  useful  in  that  the  effects  of  the 
sinusoidal  variation  can  be  directly  related  to  some  of  the  irregularities  existing  on  graphite  materials. 
More  recently,  Welsh  (9U)  and  Akin  and  Marvin  (95)  have  investigated  the  effects  of  surface  roughness 
on  spherical  geometric*  in  representative  high-pressure  environments,  Foference  95  presents  results 
obtained  in  a  combustion-driven  shock  tube.  These  results  show  an  anomalous  behavior  of  the  convective 
heat  transfer  in  that  the  data  deviate  from  the  laminar  theoretical  predictions,  increasing  directly  with 
Increasing  stagnation  presslre  (Figure  8).  The  deviation  from  theory  is  as  high  as  60  percent.  It  has 
been  postulated  that  this  deviation  is  due  to  the  flow  being  in  a  transitional  state  (from  laminar  to 
turbulent). 
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Figure  8.  Heat  Transfer  as  a  Function  of  Boundary  Layer  Thickness  (05) 
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SuBMtioo  of  the  date  presented  Id  references  80  through  95  doe*  not  provide  Day  insight 
Into  the  mechanisms  vhleh  produce  the  roughness  elements.  Knowlsdg*  of  the  materiel  processes,  groin 
die  distribution  inpurity  content ,  etc.,  coupled  vith  a  relloble  ud  detolled  theoretical  ebletlno 
model,  could  provide  this  Insight,  Id  addition,  proper  characterisation  of  the  material  constituents, 
coupled  with  a  systematic  definition  of  the  degree  r f  roughneee  required  to  produce  turbulence,  could 
result  Id  the  tailoring  Of  material*  «o  thet  1 (miner  flov  vculd  predominate  throughout  the  trajectory. 


3,3  Oranhlt*  Klcroetructure  Character! i(t ion 

Considering  the  complex  material-environmental  Interaction  poetulated  Id  (ectlon  3,3,  a  nev 
nd  dynamic  approach  to  analysing  the  ablation  performance  of  heterogeneoua  material)  li  In  progreei. 
This  approach  embodies  characterising  the  material  mlcroetructural  quantities  end  employing  the  results 
to  evaluate  individual  erosion  mechanisms  and  finally  their  combined  effect  on  high  pressure  ablation 
phenomena. 


This  eeetloo  suxata*ltes  several  of  the  aallent  features  of  this  approach  and  Illustrates  the 
methodology  employed  by  presentation  of  typical  characterisation  results .  The  illustrations  are  centered 
around  the  various  material  variables  which  may  be  of  consequence  In  the  total  ablation  behavior  of 
graphite.  The  analysis  techniques  necessary  to  ensure  a  self-conelatent  set  of  characterisation  parameters 
are  addressed  and  examples  cited  which  typify  their  application  to  Interpretation  of  the  erosion  behavior 
of  graphite.  The  data  presented  are  not  intended  to  represent  a  complete  analysis  which  le  still  In  the 
development  stage,  but  to  Indicate  the  many  facets  of  this  new  and  unique  approach  to  understanding  and 
analysing  the  ablation  behavior  of  complex  materials. 

3.3.1  Micro-  Macroporoeltv 

Porosity  of  graphitic  materials  may  Play  an  important  role  in  determining  thermochemlcal  and 
thermomechanlcal  erosion.  Considerable  effort  to  develop  reliable  methods  of  characterising  the  porcelty 
of  aerospace-grade  graphites  la  required.  This  effort  la  complicated  becauae  of  the  wide  else  range  of 
pores  that  may  occur  In  a  single  material.  In  particular,  two  types  of  porosity  are  apparent.  The  Tint 
type  arises  from  (tacking  disorder  in  the  flour  during  the  graphltlti  tion  process.  These  stacking 
disorders  produce  voids  which  may  be  only  partially  filled  by  binder  material.  The  resulting  macropores 
are  of  the  tame  order  of  site  as  the  filler  particle*.  Additionally,  mlcroporrsity  Interlaced  throughout 
the  binder  material  la  Introduced  during  graphltiaatlon  because  of  the  formation  or  shrinkage  cracks  and 
local  aeparatlon  (debonding)  of  the  binder  from  filler  partlclea.  This  porosity  is  characterised  by 
pores  in  the  aubolcron  else  range. 

An  approximate  measure  of  the  distribution  of  pore  volume  as  a  function  of  pore  radius  can  be 
obtained  from  poroalmetry  measurements.  In  these  measurements  the  graphitic  sample  is  irmerreJ  in  a 
mercury  bath  and  hydrostatic  pressure  varied.  An  effective  radius  for  the  penetrating  mercury  can  be 
determined  and  the  change  in  apparent  volume  of  the  mercury  bath  measured.  The  radius  determined  la  only 
a  qualitative  measure  of  pore  radius  since  photomicrographs  ihov  that  pores  are  generally  or  a  relatively 
large  volume  accessible  through  a  smaller  opening.  A  typical  poroelmetry  measurement  or  ATJ  graphite  is 
shown  lo  Figure  9. 

For  coarso  grain  graphitic  materials  the  accessible  volume  distributions  are  characterised  by 
two  peak*.  The  first  peak  occurs  at  very  small  pore  radius  (leas  than  0,5  micron)  end  Is  Indicative  or  the 
fln>  mlcroporea  connecting  larger  mscroporee.  The  macropore  distribution  generally  peaks  at  from  !  tu  ( 
microns  mean  pore  radius.  The  macropore  peak  can  be  virtually  eliminated  by  either  impregnation  of  the 
materiel  or  by  using  finer  filler  particles  In  the  mix. 

3.3.2  Filler  Partlclee 


The  flour  or  filler  particles  used  in  the  manufacture  of  aerospace-grade  graphite  are  important 
In  the  mechanical  erosion  model  for  graphite.  It  has  been  postulated  that  mechanical  erosion  cf  graphite 
occurs  becauae  of  direct  mechanical  removal  of  filler  particles  from  -he  oxidised  and  sublimed  surface. 
Analyses  relating  this  proceae  to  aerotherr.al  environment  have  generally  assumed  that  the  filler  particles 
can  be  characterised  by  a  alngle  mean  site  (78,  96);  however,  it  Is  more  likely  that  preferential  removal 
of  either  the  large  or  aoall  particle*  In  response  to  the  applied  pressure  and  pressure  gradient  Is  a 
strong  function  of  particle  site  and  perhaps  shape.  Hence,  the  distribution  of  particle  sizes  (Figure  in) 
may  be  a  critically  Important  factor  governing  the  mechanical  erosion  of  graphite  In  high  presaure 
environments. 


Filler  particle  site  distribution,  n(dp) 


denotes  the  number  density  of  probability  of  cccur-ence 


of  a  specific  site  filler  particle.  A  typical  curve  for  a  fine-grain  graphite  1*  presented  In  Figure  11, 
Far  tell  particular  graphite  the  distribution  is  nearly  normal.  However,  large  departures  from  Oausslan 
dist.'-b’j'ions  may  be  expected  for  coarser-graln  filler  particles.  The  mean  particle  diameter,  d  ,  and 
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variance , 


can  be  obtained  from  these  data  from  the  following  equations; 
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Additional  aharaotarl  aatl  on  of  tna  fluar  partlolaa  by  a  ahapa  factor,  A,  nay  ba  ra^ulrad  in 
aroalon  analyaaa.  Tha  ahapa  faotor  a  an  ba  dafinad  by 
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H  number  of  particles  par  gri  a 
p?  ■  filler  partial*  daneity 
a  •  parti  ole  projected  iru  par  (rut 

That*  proparti**  must  ba  measured  prior  to  th*  binder-filler  nixing  proc***.  Attempts  hav* 
bean  made  to  uaa  tha  shape  factor  aa  defined  above  In  a  relation  of  th*  fora 

A  ♦  ac  ■  t  (31) 

where  a  and  t  era  oooatant*  to  determine  th*  filler  particle  packing  density  that  could  ba  achieved  In 
a  dry  atata  (that  la,  unalxed  with  binder).  It  appear!  clear  that  no  auoh  simple  relation  exlata  to  deter¬ 
mine  th*  poroelty  of  th*  graphitic  material*  considered  for  aeroapao*  application.  In  particular,  addi¬ 
tional  parameter*  entering  Into  aa  expression  for  c  vould  Include  th*  binder  naa*  ratio,  th*  varlano*  of 
the  filler  particle  all*  distribution,  and  th*  graphitlaatlon  temperature.  Other  than  through  Influencing 
the  poroelty  of  th*  graphitic  material,  It  la  clear  at  thla  tin*  that  th*  ahap*  factor  1*  an  Important 
influence  on  erosion  neohaniam*. 

3.3.3  »indor  Properties 

The  primary  binder  properties  of  Interest  for  aorotharmal  orealon  prediction*  era  binder 
bridge  thioknaa*  and  density,  neither  Of  these  properties  can  b*  measured  directly,  hence  aealyeee  have 
been  developed  to  predict  thee*  properties  from  available  oharacterltatlon  data.  Some  ambiguity  usually 
exist*  In  the  definition  of  those  quantities, 

Th*  blndar  bridge  thickness  la  an  Important  parameter  in  determining  eroelon  rat*  due  to 
sublimation  and/or  oxidation  of  blndar  materiel  Immediately  below  t.he  eurface.  Thle  can  b*  defined  aa 
th*  mean  thioknaa*  of  blndar  material  surrounding  filler  particlaa,  macropo r*»,  and  aleropore*.  In  keeping 
with  th*  model  of  large  maoropore*  reached  by  small  openings  (mieroporea)  both  th*  filler  partial**  and 
maoropere*  era  modeled  aa  apbara*  of  uaan  volume  V  and  V  ,  respectively,  obtained  from  th*  raapactlv* 
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distribution  funotion*.  Th*  mloropora*  era  modeled  a*  cylindrical  tunnels  of  mean  radius  rBlflro  .  The 
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Inaccasalbl*  void  voIum  1*  Included  In  th*  acseeelble  aaeropor*  volume.  For  thle  model  the  mean  bridge 
thlckneee,  l.  ,  la  given  by 
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and  p  1*  the  apparent  density  of  th*  graphitic  material,  pg,ia  th*  blndar  material  density,  and  tuioro  1* 
th*  fraotlon  of  bulk  voIiom  cooslttlng  of  micropora*. 

Th*  variation  In  naan  blndar  bridge  tnicxnas*  1*  anovn  in  Figure  II  for  eev<rel  value*  of 
th*  binder  naee  fraction,  K.,  typical  of  varioui  graphite*.  Significantly,  th*  aoan  blndar  bridge 
thioknaa*  la  quit*  thlni  even  a  aaall  degradation  In  thi*  thlcknaaa  by  oxidation /sublimation  may  fra*  tb* 
otherwise  conatralnad  particle*  resulting  In  thalr  being  easily  dialodgad  from  th*  material  surfao*  by 
aarodynamlt  fore**. 

3.3. u  uma Mmiiim 

Tha  lnfluanc*  of  straaaai  on  th*  thermcmaehanioal  sroalon  can  b*  considered  from  both  th* 
macroaechanloal  and  mlcromsohanical  point  of  view.  The  result*  of  atres*  analyses  can  be  used  to  establish 
th*  aearo-  or  mloi'^failure  of  th*  material  and  relate  it  to  the  ohamlosl  and  maohanlcal  aroalon. 

In  addition  to  th*  poroalty  which  1*  praaantad  in  th*  material  Initially,  additional  surface 
and  Subaurfac*  porosity  may  ooour.  At  or  near  the  eurface  th*  Internal  microetmotur*  change*  drastically 
with  time,  temperature,  and  environment.  Correspondingly  any  change*  In  th*  materiel  properties  at  or 
near  th*  eurfao*  also  affeot  th*  macroatres***)  net  only  th*  aurfac*  atreaea*  but  aleo  tha  stress**  within 
th*  body  Itself, 

The  magnitude  of  porosity  can  be  aatlmatad  by  th*  mas*  and  heat  balance  relation.  To  a  count 
for  poroalty,  th*  material  properties  need  for  each  tlentnt  in  a  finite  element  enalyeli  can  ba  properly 
modified.  Well-Haahln  (97)  euggeatad  tha  following  modlfloatlom 
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Figure  12.  Variation  of  Binder  Bridge  Thickneei 

In  addition  to  affecting  tbe  mechanical  properties,  poroeity  ol,o  Influence,  the  egrength 
propertle*.  At  tbe  poroelty  increase*,  tbe  etrengtb  decree,#*,  Bern  experimental  reeult,  on  the  etrength 
reduction  due  to  weight  loee  (caueed  by  increased  poroiity)  of  on  AflL  graphite  are  *hovn  in  Figure  13  (98) 
Thle  etrength  reduotloo  le  predictable  theoretical!!'  by  the  following  equation  (99)  i 


o  -  1  -  f  - —  a 

p  T  e(l-K)  J 

where  ky  le  tbe  roluae  fraction  of  pore*,  IT  1*  the  volume  fraction  of  lnclutione,  (1  -  7  -  ky)  le  the 
volume  fraction  or  binder,  0  le  the  etrength  of  eolld  materiel  (ky  »  0),  and  le  the  etrength  of  poroue 
material.  8ome  typical  thaoretlcal  reaulte  on  the  effect  of  ky  and  k  on  are  ihovn  In  Figure  lk. 

Although  it  le  ioporelble  to  correlate  roeulte  given  in  Figure!  13  and  it  due  to  the  fact  that 
reference  96  did  not  give  the  Inclusion  content,  It  la  quite  obvious  that  the  theory  and  experiment  ebov 
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3  3-5  Synergistic  gff«M 

Th#  erosion  behavior  o t  graphltlo  material*  entail*  an  assessment  of  tba  oompetlng  chemical 
and  Mchanlcel  mechanisms  of  graphite  aroaloo  and  their  synerglatlc  affects.  Tba  preceding  baa  Illustrated 
eeverml  of  tba  phenomena  of  consequence  In  tty-  analytical  performance  daacrlptiont  It  raaalna  to  foraulatt 
and  experimentally  verify  tba  complete  analytical  nodal.  At  thla  stage,  It  la  claar  that  a  aora  thorough 
approach  la  necessary  la  analytical  modeling  of  hatarganaoua  materials.  Tba  benefit  la  twofold (  a 
rigorous  ablation  analyala  which  lncludaa  tha  material  alcroatructural  charactarlatlca  raaulta  and  a 
realistic  insight  to  the  aatarlal  charactarlatlca  which  contribute  to  the  observed  ablation  data  la 
obtalnad  enabling  aatarlal  formulations  to  be  made,  aimed  at  reducing  or  alleviating  tta  enhanced  erosion. 

3,U  Orsphlta  Thermo chemical  Erosion 

Theraochaalcal  erosion  of  graphitic  aatarlala,  both  In  tha  laboratory  and  In  .jrouod  teat 
faollitlea,  baa  been  the  subject  of  sane  considerable  study  for  the  past  decade.  At  present,  an  lnecaiplete 
understanding  exists  of  tha  detailed  effect*  of  material  structure  on  oxidation  and  sublimation.  It  la 
known,  however,  that  tha  surface  reactivity  of  tha  material  la  Important  In  determining  oxidation  rate* 

In  tha  Xlnetlc- coo  trolled  regime. 

Several  studies  showing  the  effect  of  graphite  solid  structure  on  oxidation  kinetics  have 
b««n  accomplished.  For  the  moat  part,  result*  are  Incomplete  because  tha  investigator*  have  failed  to 
include  all  possible  influences  on  their  experiments  in  a  controlled  mannor.  For  example,  both  the 
o.  lantatlcu  and  else  of  graphite  crystallites  have  been  neglected  in  performing  reduction  of  oxidation 
data. 


The  diffusion  cootrollad  oxidation  erosion  of  graphitic  materials  In  byperthermal  environment* 
la  relatively  well  established.  The  oxidation  kinetic  regime  has  been  formulated  by  several  investigators 
vith  various  degrees  of  sophistication.  Tha  work  by  Welsh  and  Chung  (100),  Heame  et  al.  (1C1),  sod 
Seals  (52)  provides  results  reflecting  state  of  the  art  determination  of  the  reaction  rate  regime  descrip¬ 
tion  for  graphitic  materials.  Tba  diffusion  controlled  oxidation  regime  and  transition  to  this  regime  has 
bean  described  employing  both  equal  and  unequal  cl f fusion  coefficients,  utilising  similar  boundary  layer 
assumptions  of  Lees,  the  approximations  of  Scale  and,  more  recently,  the  unequal  diffusion  coefficient 
approximation  by  Kendall  (30).  It  remains  that  significant  differences  existing  between  th*  several  acts 
of  solutions,  vhen  using  Identical  thermodynamic  properties,  have  yet  to  be  demonstrated. 

The  hlgh-tesq>*r*ture  tbermochemlcal  erosion  regime  of  graphite  Is  lee*  well  understood  because 
of  the  limitation  In  obtaining  tha  high  temperatures  required  and  associated  temperature  measurement 
errors.  Effects  of  Importance  are  vapor  pressure  of  tha  carbcn  species  evolved,  th*  possibility  of 
departure  from  phase  equilibrium.  Incomplete  definition  of  the  experiments  In  terms  of  such  considerations 
as  reaction  time*,  ate.  The  question  arises,  however,  as  to  tha  total  effect  of  tn*  idealisations  that 
era  mad#  presently  on  tn*  chemical  erosion  of  graphite.  Conclusively,  it  can  be  said  that  the  materiel 
■tructure  has  a  definite  effect  on  the  chemical  erosion  and  the  Importance  or  Its  detailed  structure 
remains  undetermined. 

3.L.1  Vapor  Properties 

Considerable  interest  and  research  have  been  focused  on  the  identification  and  thermodynamic 
characterisation  of  hlgh-temperature  carbon  vapor  species,  and  elucidation  of  the  carbon  vapor  equation 
of  state  In  the  last  several  years.  Up  until  then,  th*  thermodynamic  properties  of  all  polyatomic  carbon 
vapor  apecles  larger  than  have  bean  calculated  (102,  103,  10b)  frem  the  free  energy  Increment*  of 

Pitser  and  dementi  (102),  which  were  fitted  to  yield  a  triple-point  pressure  of  100  atm  at  LOOO°K  (105). 
Likewise,  th*  equation  of  state  calculations  (103,  10b,  106,  107,  108)  have  failed  to  include  the  heavier 
molecular  weight  species  C17  t0  C30'  despite  evidence  of  their  existence  (109,  110)  and  llkollbood  of 

their  prominence  at  high  pressure*  end  temperatures.  Th*  carbon  ring  molecule*  have  also  been  excluded. 

An  extensive  study  of  the  carbon  vapor  thermodynamics  and  thensochemlcal  ablation  performance 
of  graphite  In  hlgh-preseur*  environments  was  given  in  reference  77.  The  results  (77)  lncluled  th# 
thermodynamic  characterisation  of  blgh-tenperatur#  carbon-vapor  specie*  and  tba  elucidation  of  th*  carbon- 
vapor  equation  of  state.  Heate  of  formation  and  thermodynamic  properties  were  calculated  for  both  tba 
linear  and  mauocycllc  polygonal  ring  conformations  of  Cq  up  to  over  the  temperature  range  109°K  to 

5,000°K.  Th*  thermodynamic  properties  of  th#  linear  and  cyclic  molscule*  Cj  to  C^0  vere  calculated  by 

quantum  and  statistical  mechanical  techniques.  Th*  resonance  energy  of  each  molecule  was  computed  from 
tb*  electronic  wave  functions  by  th*  extended  Huekel  method  of  Hoffmann  (111),  and  the  corresponding 
vibrational  frequencies  were  determined  by  a  normal  coordinate  analysis  of  tho  molecular  structure.  The 
theoretical  procedures  employed  were  self-consistent  and  did  not  resort  to  th#  empirical  increments 
adopted  by  Fitter  and  Clementl  (102)  and  used  by  others  (103,  10b).  Figure  15  deplete  the  logarithms  of 
th*  partial  pressures  of  the  linear  molecules  In  phase  equilibrium  vlth  grephlt*.  The  plot  shows  linear 
Cj  is  the  most  stable  species  snd  linear  somewhat  more  (table  than  and  Cg,  In  agreement  vlth 

previous  mass  apactrooetrlc  findings  (109,  112).  It  it  also  apparent  that  the  contribution  of  ring 
nolecules  la  negligible  (77),  and  that  larger  linear  molecules  become  very  Important  at  b500°K,  and  should 
b*  included  In  calculations  performed  above  boO0eK,  In  Figure  15,  an  even-odd  alternation  is  evidenced 
above  C^0  at  U0OO®K;  the  affect  appear*  to  propagate  to  smeller  values  of  o  vlth  Increasing  temperature. 


Figure  16  reflects  the  present  state  of  affairs  In  regard  to  tola,  carbon  vapor  pressure  and 
several  calculations  that  have  been  made  In  the  past  several  years.  As  observed,  large  differences  in 
vapor  pressure  exlet  in  the  temperature  range  3000  to  b500°K.  The  calculated  total  vapor  pressure  (77) 
1*  In  good  agreement  with  th*  experimentally  determined  sublimation  point  of  Thorn  and  Wlntlov  (112)  and 
lias  between  the  recent  triple-point  temperature*  of  Schoestov  (lib)  and  Fateeva  (115).  The  calculated 
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thermodynamic  properties  given  by  Kratach  nt  al.  (T7)  ere  entirely  theoretical,  and  the  fundamental  eon* 
•tents  employed  la  the  cleculatien  were  utilised  to  establish  the  uncertainty  la  carbon  vapor  prescure  la 
a  lelf-conslstent  manner,  This  facet  of  aelf>caeelstency  ic  evaluating  the  uncertainty  In  tharmodynanls 
propartlaa  In  cerbom-apeelat  vapor  praaaura,  la  important  and  la  raflactad  by  the  uncertainty  In  tha  fraa 
anarciaa  of  formation  for  aech  apaclea  with  n  >  3. 


Tha  governing  conservation  equations  to  an  ablating  body  vtth  Internal  gaa  generation  has 
bean  atatad  In  eecticn  2.2,  Their  application  to  graphite  la  applied  In  thle  aactloe. 

Tha  elymntal  oaaa  balance  of  a  material  consisting  of  1-chaalcal  alammnta  vhoae  aurfaca 
racadaa  at  a  rata  a  and  undergoing  loaa  of  aolid  aatarlal,  mr>  at  tha  aurfaca  (Flgura  1)  la  given  hy 
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where  ■  S.  asd  a  K.  are  tha  uu  fluxes  of  tha  element  1  Into  a  control  volume  fired  at  the  receding 
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aurfaca  owing  to  an  internal  gaa  source  and  aolid  material  transfer,  respectively,  Maas  leaves  tha  control 
volume  normal  to  tha  wall  by  a  nat  flux  of  gaa,  (pv).  ,  and  hy  loaa  of  solid,  m  K,  .  Tha  mesa  flux  of  tha 

th  "  r 

1  alaaent  normal  to  tha  wall  accrues  by  a  summation  of  individual  flux  coop on ant a  consisting  of  ths 

internal  gas  sourca,  aurfaca  chemical,  reactions  with  the  external  reactive  gaa  (boundary  layer  and 
internal  gas  apaclea)  and  by  direct  aublimatlon.  For  graphite  in  an  air  boundary  layer  tha  following 
conservation  statements  era  made  at  tha  vail  (77)  t 
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Tht  gradient  of  the  individual  chemical  element*  at  the  vail  la  given  by,  e.g.,  carbon 
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and  for  unity  Lewls-Semenov  number  tha  convective  heat  transfer  to  tha  vail  by  conduction  and  diffusion 
la ,  upon  introduction  of  Stanton  number 
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Thus,  Eq,  38  can  be  aolvad  by  introducing  Eq.s  Uo,  U3,  and  bi  and,  after  rearrangement,  the 
elemental  maae  fraction  of  carbon  at  tha  wall  la  given  by 
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where  for  graphite  in  an  air  boundary  layer 
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lh«  fraction  of  aolid  entering  the  control  volume  which  la  loat  at  tba  surface  as  particulate 
solid  la  defined  from  Eq.  16 
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Inspection  of  Eq,  15  shows  that  the  mass  fraction  of  the  alemant  carbon  at  the  wall  la  deter- 
mined  solely  by  the  naas  flu*  of  elemental  carbon  Injected  at  the  surface  as  vapor,  reflected  by  the 
difference  In  fluxaa  entering  and  leaving  the  control  volvae.  Squat  loos  similar  to  Eq.  15  for  the  element 
oxygen  and  nitrogen  provide  the  necessary  relations  satisfying  Eq,  39,  A  stoichiometric  statement  then 
directly  relates  the  elemental  mass  fractions  to  the  chemical  composition  of  the  combustion-sublimation 
apeciaa  at  tha  ablating  surfaca 
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Suaerlcal  solutions  to  the  governing  equations  vara  given  In  (77)  invoking  the  chemical 
squlllbrlum  constraint.  Tha  chemical  reaction  system  included  sublimation  species  for  the  ring  and  chain 
molecules,  and  In  addition  air  apeciaa,  and  carbon-air  reaction  product*  Hated  belov 


Graphite  sublimation  ratea  over  a  large  pressure  range  are  shown  in  Figure  17  normalized  to 
the  f sailer  diffusion-controlled  oxidation  rate  (sip  ■  0.171*  Peuo  CH>.  In  the  first  several  decades  of 

pressura  tha  threshold  sublimation  tamperature  Increases  with  increasing  pressure  reflecting  the  carbon- 
vapor— pressure  diagram  shown  In  Figure  16,  The  10^  to  10^  atm  pressure  range  Is  of  special  Interest;  the 
threshold  sublimation  temparatura  Is  compressed  due  to  sublimation  of  large  quantities  of  polyatomic 
carbon-vapor  spaclss.  Thus ,  tha  high-pressure  sublimation  regime  of  graphite  occurs  over  a  fairly  narrow 
surface  temparatura  range  (circa  7000-8000°R)  In  contrast  to  Seals'*  (52)  results  (7 000-10000® F)  which 
emit  carbon  molecules  laurger  than  C  «  Upon  definition  of  e  aurfaee-ensrgy  balance  and  the  lndeptb  heat 
conduction ,  sublimation  ratea  are  directly  coupled  to  surface  temperature  end  pressure  and  the  local 
aerodynamic  beating. 

Energy  transfer  to  a  surface  undergolnlng  material  consumption  via  thermochemical  e  rot  lac  and 
thermoamchanlcal  lots  of  solid  Is  given  by  so  energy  balance.  Figure  1.  The  best  conducted  by  tbs  solid 
away  from  tha  surface  reflects  the  energy  fluxes  leaving  and  entering  the  control  volume  expressed  in  the 
following) 

q  ♦  i  .  -  a  ,  —  (pv)  H  -  m  H  ■  q  _ .  -  a  H  -  m  H  (1*9) 
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Tha  summation  in  Eq.  50  la  made  over  all  species  in  the  gas  at  th*  vail  vlth  concentrations 
defined  by  th*  mass-transfer  rstas  and  application  of  the  squlllbrlum  constant.  For  the  case  of  Le  ■  1  the 
net  heat  trsnsfar  to  tha  surface  m tr  be  axpreased  as  (at  this  tlo*  m  It  set  equal  to  zero) 
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Figure  17.  Graphite  Sublimation  Rate 
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vhen  the  usual  dafinitlon  of  the  oaa»  transfer  parameter  is  made 


Eq.  51,  when  coupled  to  a  standard  finite  difference  solution  to  the  conduction  equation  enables 
a  oanplete  description  of  the  surface  erosion  of  graphite  Including  solid  lose.  It  is  observed  thst  the 
effect  of  solid  renovel  on  the  surface  energy  balance  it  to  reduce  the  net  energy  transfer.  For  vxarple, 
at  a  surface  temperature  of  7500®?  the  B*  H  tern  le  about  1700  Btu/lb  at  100  atm  preeeure;  Eq.  51  is 

reduced  only  by  850  Btu/lb  for  f(  «  0.5.  In  comparison,  the  enthalpy  of  the  combustion-sublimation 

products  at  the  ablating  surface  le  about  «00  Btu/lb  as  shown  in  Figure  18.  the  enthalpy  of  the  gaseous 
products  (given  by  Eq.  50}  evidences  ths  sxothsralc  CO  chemical  rsactlon  product  belov  5000°?  surface 
temperature.  The  increase  In  enthalpy  abova  5000®?  is  coincident  with  aubllmatlon;  highly  endothermic 
carbon  vapor  species  svolvs  at  ths  surfsce. 

Details!  results  of  ths  high-pressure  sublimation  regime  are  shown  in  Figure  19.  The  chemical 
composition  at  ths  ablstlng  graphite  surfsce  is  given  st  100  atm  pressure.  In  addition  to  the  air-carbon 
reaction  products,  a  complex  set  of  sublimation  species  is  evidenced,  Ths  triatoaie  carbon-molecule  is 
most  ab undent ;  however,  s  lerge  number  of  polyatomic  molecules  of  rather  low  concentration  are  seen. 

Other  them  C-  the  prevalence  of  any  single  carbon  vapor  spades  Is  not  evident,  Although  vapor  species 
baring  ring  structures  wars  included  in  ths  numerical  computations,  the  concentrations  shown  are  due  only 
to  linear  structured  molecules.  Conclusively,  ths  ring  molecules  can  be  safely  neglected  In  performance 
analyses  for  graphitic  materials.  Omission  of  the  linear  carbon  species  at  high  pressures  examined  hers 
may  be  quits  critical  as  affects  mass  transfer  (77). 

It  Is  interesting  to  examine  how  ths  ablated  mass  of  carbon  la  partitioned  among  the  chemical 
composition  shown  in  Figure  19.  In  Figure  30  the  elemental  mass  fraction  of  carbon  is  shown  extending  from 
ths  diffusion-controlled  oxideton  regiue  well  into  the  sublimation  regime.  As  indicated,  the  constant  mats 
fraction  of  0,15  raflscta  ths  carbon  mast  in  CO.  At  about  4500°?  nitrogen  reads  forming  tbs  cyano  radical 


Jk 

•ad  cyanogen  i  only  *  smell  of  carbon  is  ablated  because  of  these  reaction  products.  At  about  tbs 
ssm  temperature  (6500®F)  contributions  dus  to  C;  to  Cj  occur,  increasing  ths  mess  fraction  of  eblnted 

carbon.  With  increased  temperature  larger  carbon  vapor  spselss  sublime  from  ths  surfsco  appearing  first 
ss  spool »s  having  fros  It  to  10  carbon  atoms  sod  than  In  spselss  h trine  fras  11  to  20  carbon  atoms.  For 
ths  100-eta  prsssurs  condition  shown ,  llttls  contribution  from  vapor  apscl  ss  with  greater  than  20  carbon 
atoa*  occurs.  Mors  clearly  than  indicatsd  in  Figure  19 ,  it  is  apparent  that  ths  additive  snail  concen¬ 
trations  of  ths  polyatomic  carbon  vapor  spselss  la  quits  appreciable. 


Figure  20.  Partition  of  Elemental  Mass  Fraction  of  Carbon  with  Surface  Temperature 


3.*t.3  Sensitivity  Analysis 

The  uncart airty  limits  established  for  ths  carbon-vapor  species  as  affacts  mass  transfer  and 
energy  transfer  can  be  assessed  employing  the  conservation  aquations  previously  described.  Figure  21 
tbs  normalised  sublimation  rata  for  several  pressures  included  *20  percent  uncertainty  in  fras  eesrjy  of 
formation  for  ths  carbon-vapor  spselss.  The  uncertainty  band  rsflactsd  in  Figure  21  is  consistent  with 
the  vapor  pressure  diagram  shown  in  Figure  16^  the  higher  pressure  evidences  a  greater  uncertainty  band 
dus  to  the  greater  concentration  of  aubllmatlou  specias.  For  comparative  purpose!,  results  calculated  by 
3eala  (52)  are  shewn  for  a  normalised  sublimation  rate  of  li.O.  Also  shown  are  results  by  Dolton  et  el. 
vlOb)  at  1  aim  pressure.  He  includee  vapor  specias  up  to  C^,  employing  empirically  determined  free 

energies  forced  to  fit  a  presumed  Uo00°K  triple-point  temperature.  It  is  teen  that,  at  100  atm  pressure 
and  for  cenatant  surface  temperature,  sublimation  rats  can  be  as  much  as  3  times  larger  than  ths  nominal 
value.  TTie  effect  of  uncertainty  in  thermodynamic  propart  Isa,  diminishes  with  decreased  pressure,  owing 
to  ths  smaller  vapor-pressure  contribution  of  the  high  molecular  weight  vapor  products.  Corresponding 
results  shoving  ths  effect  of  ths  combustion-sublimation  products  an  wall  enthalpy,  are  shown  in  Figure  22. 
The  coaposlt*  of  these  results  enable  an  assessment  of  ths  uncertainty  on  surface  recession  rate  to  be 
made. 


The  calculated  surface  recession  rate  for  e  range  of  environmental  conditions  studied  in 
correlation  of  high-pressure  teet  date  (77)  is  shown  in  Figure  23.  At  the  highest  stagnation  enthalpy 
shown  (10,000  Btu/lb)  calculated  surface  temperature  is  T500t>B  at  100  atm  prsssure.  Ths  sensitivity  of 
surface  receeslon,  owing  to  the  uncertainty  band  on  properties,  is  snail.  This  sensitivity  decreases  with 
both  reduction  in  stegnetlon  enthalpy  and  pressure,  evidencing  the  effect  shown  in  Figure  21  cm  sublimation 
rate.  Even  though  s  large  uncertainty  in  the  thermodynamic  properties  of  the  carbon  vapor  species  was 
included  (resulting  in  almost  a  decade  uncertainty  lu  carbon-vapor  pressure)  the  effect  on  surface  recession 
is  less  than  10  percent  within  the  environmental  conditions  diseuassd.  This  result  occurs  bscsute  of  the 
counter-balancing  affect  of  increased  endothermic  heat  absorption  when  greater  quentlttes  of  carbon  vapor 
are  evolved  at  tbe  surface.  Thus,  when  using  s  sslf-cooslstent  model  for  calculating  the  vapor  species 
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Figure  23.  Themochemicai  Surface  Secession  Pete 


thermodynamic  properties  with  eseocieted  uncert*int>  .  .ivlty  o  surface  recession  is  comparable 

to  the  accuracy  for  which  local  aerodynamic  heating  cu.  ions  can  be  define-.. 


alts  Mechanical  Erosion 


This  section  presents  the  mechanical  troalon  aapecta  of  graphite  given  in  (77)  and  it  devoted 
to  determination  of  the  fraction  of  particulate  solid  eroded  from  the  receding  surface  in  terms  of 
environmental  exposure  condition*.  To  accomplish  identification  of  the  thsracaecbanlcal  recession  com¬ 
ponent  from  aurfae*  erotic*  m**sure«iants  of  the  test  data,  the  thsrmocbemlcal  recession  data  described  in 
section  3.*  i*  used.  The  probable  wseertainty  in  thermodynamic  data  employed  on  turf see  rvceeaiou  rate*  was 
evaluated  in  a  description  of  aublipatlon  rates.  Lots  of  partieulato  solid  st  the  heated  surfec*  is 
addressed  in  this  section  as  affects  mass  and  energy  transfer.  Correistion  of  the  experimental  date  in 
terms  of  thermocbamlcal  erosion  enables  quantitative  identification  of  the  solid  fraction  eroded. 


3.5.1  Ablation  Data 


Representative  graphite  plasma  arc  ablation  data  given  in  (77)  lo  presented  for  completeneee. 
A  typical  block  of  date  is  shown  in  Figure  2lj  experimental  recession  date  are  plotted  against  apparent 
stagnation  enthalpy.  Two  set*  of  data  ere  given  as  typical.  Shown  are  data  taken  in  the  etagnatlon 
pressure  range  of  83.5  to  68.6  eta  end  date  gathered  between  18.6  end  23.1  etm.  Eecb  set  of  data  was 
correlated  (77)  using  e  least  square »  curve  fit  of  the  data  point*.  Results  similar  to  the**  have  beeo 
obtained  fur  ncetlnal  stagnation  pressures  of  55  end  70  etm.  The  composite  of  the  recession  date  et.eclee 
construction  of  the  experimental  recession  rate  map  ehcrwn  in  Figure  25.  Here  experimental  recession  rat* 
i*  shown  in  terms  of  stagnation  pressure  end  stagnation  enthalpy.  Figure  25  presents  the  total  picture 
of  the  dependency  of  the  experimental  receesiou  rets  data  upen  the  environmental  condition*,  stagnation 
enthalpy,  and  pressure .  The  data  appear  to  be  linear  function*  of  pressure  end  enthalpy. 


Ccnstructiu..  of  the  experimental  recession  data  map  enable*  a  comparison  of  the  calculated 
thsrmocbemlcal  recession  to  the  obaerved  experlmantal  raeeesion  data,  rigur*  6  presented  the  ratio  of 
experimental  and  thermochemical  mas*  fluxes  for  the  teet  environmental  spectrum  shown  in  Figure  25. 


Figure  »6  ahowa  tba  daparturt  tram  tbaraoabtaiaal  rroalon  In  aort  aaanlngful  tint,  from 
tba  au^siita  of  aaparlaantal  data  tad  nuaarloal  raaulta,  lueluding  aolld  raaowal,  tha  fraction  of  *olld 
arodad  ka«  km  aalaulatad  aad  ii  daplatad  la  Flgur#  26.  Tha  tbraahold  atafnatire  praaaura  at  vhleb  lota 

of  parti  aula  to  aeUd  iaitiataa  la  daposdaat  on  atagaatioo  aatbalpjr  (raflaetioa  aanvaatlva  boat  rata  at 

aaaataat  praaaura).  lha  raaulta  abown  art  oonalatrnt  with  thoaa  arpaotad,  laarataad  at  am  at  ion  enthalpy 
raaulta  la  hlfbar  aurfaaa  taaporaturaa (  henna,  greater  aubturftat  loaa  of  bladar  aattrlal  radueio* 
k  India  |  af  t  lUor  part  la  la  a  within  tha  graphlta  aatrlx.  Xoeraaainf  praaaura  above  tba  tbraabold  value 
aauaaa  a  large  laaraaaa  la  tba  aolld  fraatlon  eroded  orar  a  fairly  a  arrow  praaaura  ran  (a.  Tha  loarataant 

la  aalld  arodad  tbaraafta.  daaraaaaa  aa  tba  full  praaaura  rang*  la  aovarad. 


flgur*  26,  Fraction  of  dolll  Krodad  Fro*  Cnparioaatal  Data 

Caaplanentary  raaulta  art  (Ivan  in  Flgur*  2T  dapietini  tba  pradlatad  aurfaaa  taaparaturt  aa 
affaatad  by  laaa  of  aalld  from  tha  haatad  aurfaea.  lolavorthy,  ia  tha  rathar  aaall  daaraaaa  in  aurfaaa 
tenperttur*  for  avan  laria  fraatlona  of  aolld  arodad.  Snrfao#  taaparatura  la  aoat  affaotad  by  aolld 
lata  at  tha  low  et agnation  enthalpy i  tha  hl|h  atagnatlon  anthalpy  aomvhat  aaaka  tbo  laaa  offlolaot 
romval  of  Mtorial  in  tbo  aolld  atata, 

To  lntarprat  tha  raaulta  In  Flcura  26  In  tame  of  tha  oontanta  of  tbia  ravlav,  rafaranoa  la 
and*  to  tha  poatulatod  prafartntlal  aroalon  aaohanltaa  provJoualy  dtaouauad.  Flrat,  it  vaa  obaarvod  that 
tha  affaet  of  praaaura  on  aurfaoa  toaparaturo  la  aaall  at  tbo  raapaotlva  valuta  of  atafnatlon  anthalpy 
(within  aha  praaaura  rant#  aovarad).  Thua,  oath  valut  of  atagnatlon  anthalpy  aan  bo  viawad  lu  taraa  of 
naarly  oonatuit  aurfaea  '•■parature  for  tha  rani#  of  aolld  fraatlon  arodad  aridancad  in  Fiiuro  26. 

Conaldtr  than  tba  rola  of  « tarnation  praaaura,  tha  nornal  fare*  aatlni  on  tha  haatad  aurfaea.  Aa  tba 
aurfaaa  raaadaa,  flllar  partial#!  of  varloua  alia  ara  oontlnuoualy  axpoead.  Loaa  of  binder  balov  tha 
haatad  aurfaaa  and  by  direct  tharaoahaalaal  aroalon  at  t)«  aurfaaa  tanda  to  fraa  tha  otharvlat  aonatralnad 
and  aoMvbat  Interleaved  parttolaa.  Tha  largar  flllar  partial#*  ara  than  rathar  aaaily  blown  fraa  by 
•v*a  aodarat*  praaaura*,  eomanelne  vlth  tha  hlihor  traporatur*  axpoauxa  jndltlot,*  (I. a.,  atafnatlon 
anthalpy)  with  lnarttaad  praaaura  aatlni  to  fraa  tha  aaallar  partlalaa.  Alio,  alnuta  aaounta  of  blndar 
laava  tha  aurfaea  attaohad  to  tha  dlalodfod  flllar  partiolaa.  Aa  Indicated  In  Fl|ura  26  tha  affaot  of 
anthalpy  (hanoa,  aurfaaa  tmparatura)  at  tha  hl|har  praaaura*  appear*  to  dlalnlah  and  praaaura  doalnata* 
tha  loaa  of  pai-tlaulata  aolld,  Intareatlnily ,  tha  fraatlon  of  aolld  arodad  appear*  to  b*  approMhloi  an 
aayuptotls  ralu*  in  tha  rani*  0.6  to  0.6,  liinlflaantly ,  tba  noaloal  m*i  fraatlon  of  filler  partiolaa 
(about  0.79)  fall*  within  thia  rani*.  Finally,  tb*  atfao*  of  aarod/naalo  ahaar  and  praaaura  (radiant* 
typtaal  of  pofltion*  roaovad  from  tha  (tarnation  point  la  to  tnduo*  inaraaaad  partial*  loaa.  Thua,  on* 


Figure  27.  Influeno#  of  Solid  Fraction  Iroded  on  Burrao*  Temperature 

la  fanad  with  loaa  of  particulate  graphite  ovar  a  ratbar  broad  rang#  of  aerodynamic  conditions,  par hap a 
avail  extending  into  tha  relatively  moderate  euper-orbltel  antiy  anvlromant. 

3.6  £flBai.uj»Vga.t.Hfla  £UMmfM\9R» 

A  nodal  incorporating  tha  allotroplo  faaturae  of  bulk  graphite  bae  been  poetulated  and  used  to 
describe  tbo  continuous  removal  of  particulate  graphite  from  the  heated  surface.  Although  hypothetical  in 
its  formation,  sufficient  laboratory  data  ars  advanced  to  establish  tba  creditability  of  removal  of  filler 
particles  in  tba  aolid  state.  Tha  key  role  of  the  different  forms  of  carbon  in  graphite  on  nurfaoe  erosion 
has  bean  established.  Thsraodynaalo  properties  of  the  linear  and  cyclic  molecules  C3  to  C30  vers  calculated 
by  quantum  and  statistical  mechanical  techniques.  The  resonance  energy  of  each  molecule  vas  computed  from 
the  respective  electronic  vara  functions,  and  the  corresponding  vibrutic.i  frequencies  vers  detnrmioed  by  e 
normal  coordlnatt  analysis  of  tha  aclscular  structure.  The  theoretical  procedures  are  eelf-coneietent , 
retort  to  tbs  empirical  Increments  sdoptsd  by  Fitter  and  Cleaenti  and  used  by  others  has  been  avoided.  The 
agreement  of  the  oalouleted  vapor  pressure  with  the  triple-point  and  sublimation  data  la  excellent.  By 
purely  analytical  naans,  tha  uncertainty  in  tha  fras  enargy  of  formation  has  it*n  established;  the  effect 
cm  oaloulated  thermochemical  surface  recession  vas  found  to  be  small.  Significantly,  tbls  result  has  en¬ 
abled  identification  of  the  solid  fraction  of  graphite  eroded  in  high-pressure  ground  taste  vith  confidence. 

To  establish  the  solid  fraction  of  particulate  graphite  eroded  ablation  data  encompassing 
fine  to  coarse  grain  graphites  have  been  correlated.  The  loss  of  materiel  in  tbs  solid  stats  bss  been 
determined  over  a  range  of  environmental  conditions  (in  terms  of  stagnation  pressure  and  enthalpy) 
evidencing  large  fractions  of  solid  removal.  ix>ss  of  particulate  solid  Initiates  st  rstbsr  lov  stagnation 
pressures  1  the  actual  threshold  is  dependent  on  hasting  rate  as  veil  as  pressure.  The  mechanisms  by 
vhioh  lose  of  partloulats  graphite  oocurs  appears  to  ba  by  preferential  surface  recession  and  subsurface 
loss  of  binder  matrix.  Filler  particles  ere  then  easily  dislodged  from  the  surface  by  even  moderate 
foroes  (l.s.,  pressure,  sheer).  Surprisingly,  the  loss  of  solid  st  the  heated  surface  doss  not  appreciably 
depress  surface  temperatures,  even  for  rather  large  fractions  of  solid  loss. 

The  results  prassntsd  in  this  rsvlsv  (77)  open  several  interesting  avenues  for  pursuit. 

Clearly,  the  role  of  graphite  microstructura  In  dictating  surface  erosion  cannot  be  dismissed.  Closely 
related  to  the  mlcrostruoturs  of  tl.e  binder  and  filler  phases  ars  the  thermophyelcol ,  mechanical,  and 
obemleal  properties  which  individually,  or  In  combination,  affect  ablation  performance.  Future  analyses 
employing  statletioal  distributions  of  particle  sics,  void  sits  and  spacing,  and  accounting  for  the 
temperature-time  oouree  of  the  crystallite  growth  will  be  necessary  before  a  completely  analytical  per¬ 
formance  description  of  graphite  can  be  achieved.  Finally,  nev  graphites  overcoming  the  evident  dissimilar 
characteristics  of  the  constituent  material*  appear  to  afford  a  means  of  reducing  or  alleviating 
preferential  erosion,  hence,  lose  of  solid  material. 
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4.0  HATtRlAtC  BEHAVIOR  ABO  TESTING 

The  Importance  of  materials  character  it  at  Ion  In  tha  study  of  ablation  phenomena  cannot  be  over 
emphasized.  Slnca  Boat  appllcatlona  require  specialized  naterlala  and  alnca  each  envlronBent  lotaracta 
differently  vlth  the  ablator,  It  may  be  aaid  that  each  ablation  theory  la  fathered  by  development  of  an 
lap roved  ablator.  It  la  dealrable  to  revlev  factor*  which  affect  the  ablation  performance  of  material, 
eeneltlvlty  of  ablation  to  varloua  material  and  environment  parameter!  and  maana  of  ablation  teating. 


The 

materiel  classes 

of  Interest  are: 

1. 

Refractories  and 

ceramics 

2. 

Plastics 

3. 

Elastomer* 

4. 

Composites. 

Energy  of  deeompoeltlon,  bum*  tranafer  blocking  of  convective  heating,  thermal  dlffuaivity, 
and  aurface  ablation  temperature  are  acme  of  the  ablation  parameters  Influencing  the  ablative  material 
aelection.  Table  3  summarizes  application  of  ablative  material!  and  tha  general  Importance  of  selection 
criteria.  The  aurface  ablation  tsmperacures  are  Indicated  to  be  high,  moderate,  or  lov  at  characteristic 
of  efficient  ablative  materials  for  the  application.  Tha  heating  exposure  time  strongly  Influences 
lnsuletlve  requirement! . 

The  follovlng  paragraphs  Identify  the  sensitivity  of  ablation  parameters  to  material  properties, 
environmental  effects,  and  thermo* tructur el  consideration*. 

4.1  Material  Property  Effects 

Material  properties  have  a  broad  influence  on  mass  ablation  (and  surface  recession)  and  on 
thermal  loeulatlve  efficiency.  The  thermal  insulation  requirements  era  defined  by  the  temperature  which 
the  structure  and  heat  shield  bond  are  capable  of  vlthatandlng.  Total  weight  requirements  result  from 
the  anticipated  mas*  lost  and  insulation  demand. 

The  maas  ablation  comparisons  In  Figure  28  show  that ,  although  a  carbonaceous  material  la 
aor*  susceptible  to  oxidation,  total  mass  lose  la  greater  for  the  elllca  materials.  The  major  reason  la 
that  aurface  ablation  temperature*  are  greater  for  carbons,  thus  Increasing  surface  reradlatlon  and 
reducing  the  convective  heating  -  also  t  «  0.9  for  carbon  base  materials  and  e  «  0.5  for  silica  base 
materials.  8 
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TABLE 


ABLATIVE  MATERIAL  APPLICATIONS  AID  SELECT! OH  CRITERIA 


Application 

Heating 

Exposure 

Time 

(Seconds) 

Weight  Per 
Unit  Surface 
Area 

Maas 

Ablation 

r>ata 

— 

Thermal 

Insulatlve 

Properties 

Surface 

Ablation 

Temperature 

Surface 

Recession 

1. 

Hlgb-Parforaanca 
Conical  Raantry 
Vahiclas 
(v/cnA  >  1500 

D  P»f): 

15  to  10 

a.  Hoaacap 

Unimportant 

Um 

Important 

Unimportant 

High 

Important 

b.  Conical 

Sactlon 

Important 

Important 

Important 

High 

Use 

Important 

o.  Aft  Covor 

Important 

Unimportant 

Important 

Low  to 
Moderate 

Unimportant 

2. 

Lov-Balllatlc 
Cosff lol.nt 
Raantry  Vanlclat 
(w/c„A  <  1500 

V  P»f): 

20  to  100 

a.  Hoaacap 

Lesa 

Important 

Important 

Less 

Important 

High 

Important 

b.  Other  Raglona 

Important 

Important 

Important 

Moderate 
to  High 

Unimportant 

e.  Aft  Coaar 

Important 

Unimportant 

Important 

Lov 

Unimportant 

3. 

Low  to  Inter - 

medlat#  Range 

Ballistic 

Ml*. lie. 

10  to  50 

a.  Hoaacap 

Lass 

Important 

Important 

Lesa 

Important 

Moderate 
to  High 

Lesa 

Important 

b.  Otner  Raglona 

Important 

Important 

Important 

Low  to 
Moderate 

Unimportant 

c.  Aft  Cover 

Important 

Unimportant 

Important 

Low 

Unimportant 

U. 

Lifting, 

Satellite  and 
Planetary  Entry 

300  to 

2000 

Important 

Unimportant 

Important 

Moderate 

Important 

5. 

Launch  and 

Booster 

Insulation 

5  to  200 

Important 

Unimportant 

Important 

Low 

Unimportant 

o. 

Rocket  Nozzles 

1  to  100 

a.  Throat 

Less 

Important 

Unimportant 

Lets 

Important 

High 

Important 

b.  Other  Regions 

Important 

Less 

Important 

Important 

High 

Lass 

Important 

7. 

Hypersonic  Inter¬ 
ceptor  Missiles 

5  to  100 

a.  Hosecap  and 
Control  Pin 
Leading  Edges 

Unimportant 

Leas 

Important 

Leas 

Important 

High 

Important 

b.  Control  Pins 
Forward  Sur¬ 
face  Regions 

Important 

Important 

Important 

Moderate 
to  High 

Unimportant 

c.  Other  Regions 

Important 

Le.s 

Important 

Important 

Low  to 
Moderetn 

Unimportant 

figure  28 


Km*  ablation  rate*  era  enhance!  by  llqull  flov  ablation,  ipallation  of  the  char  layer,  and 
mechanical  removal  of  the  char  layer.  Silica  base  material*  are  eueceptlble  to  liquid  flov  eurface 
receeelon.  However,  It  la  only  important  for  nosecap  application*  and  other  region*  subjected  to  elgnl- 
flcant  pressure  and  shear  gradients.  A  detrimental  characteristic  of  many  charring  ablators  1*  lov-char 
strength.  The  application  of  ablative  materials,  such  as  the  epoxies  and  elastomeric* ,  is  consequently 
restricted  to  low-shear  and  lav-pressure  regions  or  to  vehicles  subjected  to  less  severe  heating  environ¬ 
ments.  Shear  or  spallation  of  the  char  is  also  an  important  consideration  to  high  strength  charring 
ablatora . 

Maaa  ablation  may  not  ba  critical  for  some  applications.  Ablative  material  selection  la 
strongly  sensitive  to  material  properties  which  Influence  thermal  inaulative  performance.  A  numerical 
dsicriptios  of  the  transient  heat  conduction  relationship  indicates  the  sensitivity  to  aatsrial  properties 
with  the  dominant  material  property  being  the  thermel  conductivity  since  the  specific  hsete  for  most 
ablative  materials  hsve  nearly  identical  values. 


Heat  coaduetad  into  the  ablative  material  it  absorbed  by  pyrolysia  decomposition  and  by  the 
naae  transfer  toward  the  eurface.  The  decreased  density  due  to  internal  pyrolysis  results  in  the  thermel 
conductivity  shifting  from  low-undecomposed  material  value*  to  higher  values  for  the  char.  An  example  of 
the  influence  of  high-char  thermal  conductivity  on  phenolic-nylon  ablation  ie  shovn  in  figure  29.  The 
thermal  conductivity  data  are  from  (117).  Although  surface  recession  is  lov,  the  buildup  of  the  char 
thlckaeaa  near  the  end  of  flight  1«  about  0.5  Inch.  The  char  denaity  la  about  0.30  tinea  the  original 
density  of  72  lb/ftJ;  consequently.  Internal  mass  transfer  removed  a  significant  fraction  of  heat  con¬ 
ducted  in  at  the  surface. 
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Figure  29.  Surface-Recession,  Pyrolysis-Depth,  ud  Surface-Temperature 
Histories  for  Phenolic  nylon  Ablation 

The  lnsulatlve  performance  of  ablation  material*  la  a  atrong  function  of  thermal  conductivity 
aa  seen  in  the  above  example  and  (118).  The  thermal  conductivities  of  several  ablation  materials  are 
presented  In  Figure  23  as  functions  of  temperature.  The  experimental  correlation  for  phenolic  nylon  (117) 
is  compared  vlth  an  unpublished  correlation  derived  with  a  charring  ablator  computer  program  vhlch  van 
used  to  correlate  experimental  ablation  data.  The  data  shovn  In  Figure  30  are  representative  of  different 
classes  of  ablative  materials.  Tungsten  represecta  the  refractory  metals^  graphites  represent  a  class  of 
materials  by  themselves,  boron  nitride  is  a  representative  refractory  ceramic,  carbon  phenolic  le  a  high- 
temperature  ana  char  strength  composite,  and  phenolic  nylon  1s  representative  of  the  plaatlc  composites. 
Purple  blend  MOD- 5  Is  one  of  a  class  of  materials  specifically  tailored  to  provide  lov-thermal  conductivity. 
Hydrogen  gas  and  air  thermal  conductivities  are  shovn  for  comparison  purposes. 

For  applications  requiring  good  lnsul&tlve  properties  coupled  to  long  heating  exposure  time, 
namely,  lifting  reentry  vehicles,  ablative  plaatlc  Or  elastomeric  composites  have  been  developed  vlth  the 
follovlng  properties  to  give  improved  lnsulatlve  and  ablative  performance:  (a)  surface  oxidation  resistance, 
(b)  lav  density,  (c)  lov-ebar  density,  (d)  lov-oolecular  weight  pyrolysis  products,  (e)  fiber  reinforcement 
for  char  strength,  and  (f)  spherical  voids.  Examples  of  these  Materials  are  purple  bleDd  MOD-;  and 
phenolic  cork. 

Pyrolysis  decomposition  predicted  using  the  Arrhenius  rate  equation  has  limited  effects  on 
ablative  material  performance  provided  the  decomposition  temperature  Is  adequately  duplicated,  l.e.,  from 
about  800°R  to  150C°R. 

For  materials  vhlch  form  thick  char*  evidence  of  "cracking"  of  gaseous  hydrocarbons  hat  been 
found  through  observed  increases  In  densities  of  chars  near  the  ablating  surface.  The  effect  on  the 
ablative  material  char  le  to  Increase  Its  strength  and  at  the  same  time  to  increase  the  char  thermal 
conductivity.  Char  layer  spallation  can  result  from  deposition  causing  decreases  In  char  permeability 
and  increases  In  internal  pressure. 

Environmental  Effects 

As  Indicated  above,  ablation  performance  is  sensitive  to  environmental  effects.  An  Increase 
In  the  convective  heating  rate  usually  Increases  ablation  rates.  However ,  If  there  Is,  e.g.,  a  proportion¬ 
ate  increase  In  reradlatlve  effects,  the  Increase  vlll  not  materialise.  An  example  of  this  effect  Is  an 
Increase  In  heating  caused  by  higher  enthalpy.  Such  changes  In  beating  Increase  the  ablation  efficiency 
of  the  materials  (ll8),  and  therefore  do  not  Increase  ablation  rate*  observed. 

One  effect  of  boundary  layer  enthalpy  potential  on  ablation  efficiency  is  in  ccsibuation 
reactions,  vhlch  can  be  accounted  for  (119)  by 


In  the  follovlng  relationship 
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Figure  3C.  Comparisons  of  Thsrmal  Conductivltiea 

where 


»  heat  in  combustion  per  unit  mats  of  oxygen 
K^#  »  oa«e  fraction  of  oxygen  at  the  outer  edge  of  the  boundary  layer 

The  term  q^  gives  an  upper  hound  on  the  effect  of  combustion  by  assuming  the  all  oxygen 

diffusing  toward  the  surface  reacts  with  the  material  vapor.  For  silica  base  materials  this  is  not  true 
and  in  a  rocket  nozzle  environmental  oxygen  may  be  fixed  as  carbon  monoxide. 

In  Eq  (53)  the  combustion  contribution  (Koeh  )  to  heating  and  to  the  mass  transfer  blocking  la 

independent  of  the  boundary  layer  enthelpy  potential.  Consequently,  combustion  heating  reactions  are 
proportionately  less  Important  for  high  enthalpy  potectlal  heating  environments  that  reentry  vehicles 
experience  over  major  portion*  of  their  flight.  The  variation  in  material  efficiency  that  results  with 
an  Increased  enthalpy  is  shown  In  (113). 

Environmental  pressure  and  shear  forces  strongly  influence  ablation  performance  for  reentry 
vehicles  and  to  a  lesser  extent  for  throats  of  rocket  nozzles.  The  charring  ablators  are  recognized  to 
have  char  strengths  considerably  lower  than  the  strength  of  the  undecompoeed  materials.  Consequently,  those 
materials  which  form  weak  chars  are  restricted  to  application  in  low  pressure  and  shear  environments. 


It.  3  Tettlax 


Ablation  tsatlng  It  the  eethoc"  by  which  the  analytical  model!  are  verified,  many  material 
properties  are  deterained,  and  candidate  heat  shield  materials  are  screened  and  selected.  Because  of  the 
importance  of  testing  in  providing  preflight  design  verification  and  confidence,  it  is  obviously  desirable 
to  perform  these  tests  lu  an  accurately  simulated  if  not  duplicated  flight  environment.  For  many  reasons 
(primarily  facility  power  requirements  and  model  site  limitations)  this  la  usually  impossible,  and  tests 
smst  be  run  under  several  types  of  partial  flight  simulation.  These  results  can  then  be  evaluated  by 
mean*  of  pravlously  proven  mathematical  ablation  models,  or  by  overlapping  the  teets  in  such  a  way  that  a 
coherent  composite  picture  of  the  important  phenomena  can  be  constructed  (12). 

Ballistic  and  lifting  entry  vehicles,  designed  for  various  missions,  experience  a  vide  range 
of  entry  environments.  Two  convenient  paremetera  for  Illustrating  these  environments,  vhlch  depend  on 
flight  performance  rather  than  on  geometry,  are  total  anthalpy  and  stagnation  point  pressure.  Figure  31 
is  e  plot  of  pressure-enthalpy  with  overlays  of  typical  entry-vehicle  trajectories  and  maps  of  approximate 
simulation  facility  performance.  Two  important  conclusions  may  be  immediately  reached: 

1.  Mo  one  facility  can  duplicate  a  complete  flight  environment,  although  large  portions  of 
some  lifting-body  trajectories  can  be  closely  approximated. 


2.  High  pressure  and  high-enthalpy  effects  cannot  be  duplicated  simultaneously. 


Figure  31.  Comparison  of  Flight  and  Ablation  Teat  Facility  Performance  Parameters 

For  simulation,  steady  state  ablation  is  usually  desired  since  this  condition  provides  tbe 
necessary  materiel  parameters,  evaluation  for  tbe  analytical  models.  For  this  reason  the  facilities  most 
applicable  to  entry  eblatlon  studies  are  the  plasma-arc  tunnel  and  the  rocket-motor  exhaust,  both  of  vhlch 
have  operating  times  ranging  from  several  seconds  to  several  minutes.  Facilities  such  as  the  ballistic 
range  and  shock  tube  are  of  only  slight  interest  because  of  their  very  short  test  times  —  approximately 
10  milliseconds  and  100  microseconds,  respectively. 

Because  of  simulation  difficulties,  test  conditions  oust  be  carefully  selected  to  obtain 
meaningful  results.  Plasma-are  facilities  can  frequently  provide  adequate  simulation  up  to  entry  velocities 
of  30,000  ft/sec.  The  important  simulation  parameters  for  most  ablation  performance  and  screening  teats 
are  the  beat  flux  and  tbe  enthalpy  (121). 


Efficient  performance  of  the  glassy  ablators  depends  on  the  retention  of  the  liiuld  melt  layer, 
to  tht  aerodynamic  shear  levels  of  flight  should  also  be  duplicated,  if  possible.  The  ablation  performance 
of  charring  ablators  is  known  to  be  reduced  under  very  high  pressures  or  pressure-gradients  vhicb  tend  to 
mechanically  remove  the  char  layer.  For  flight  vehicles  which  will  experience  unusually  high  local 
pressures,  such  as  low-drag  ballistic  entry  vehicles,  it  may  be  necessary  to  simulate  the  pressure  and 
shear  levals  as  veil  as  tbe  beat  flux  and  enthalpy.  However,  since  this  amounts  to  full  scale  flight 
duplication,  these  condition!  must  be  simulated  in  separate  tests  and  applied  with  much  caution  tc  the 
flight  environments. 


1.6 

At  the  very  high  «ntry  velocities  which  will  be  experienced  In  interplanetary  return  mlselons, 
the  redlant  beating:  free  the  ihock-he eted  air  een  be  equal  to  or  significantly  greeter  then  the  convective 
heetlng  retee.  Since  the  efficiency  of  most  ebletlve  materials  depend*  on  the  heet  blockage  effect  of  the 
transpiring  genes  which  will  have  little  or  no  effect  on  the  redient  besting,  ebletlon  performance  will  be 
reduced.  Boms  plntae-ere  fscllltles  ere  erelleble  (122)  which  cen  superimpose  redient  heetlng  on  the 
convective  heetlng,  end  these  must  be  employed  for  performance  evaluation  of  materials  to  he  uted  In  this 
type  of  environment. 

Because  of  the  difficulties  sxperlenced  In  simulation  of  flight  conditions  in  plasma  ere 
facilities  end  variation*  In  tupposedly  similar  test  conditions  from  facility  to  facility,  several  "round 
robin"  test  programs  (123,  12b ,  12$,  126)  were  conducted  to  compere  the  facilities  end  allow  for  a  sore 
meaningful  comparison  of  test  results.  Also,  several  comparative  materiel  performance  tests  (126,  12T) 
for  seme  materials  were  conducted  In  a  "round  robin"  fashion.  As  a  result,  better  comparisons  are 
available.  Characteristics  of  various  facilities  are  also  given  in  the  above  references,  as  well  as  in 
(128,  129),  and  vill  not  be  tabulated  again  in  this  review.  Due  to  the  periodic  changes  Id  existing 
facilities  range  of  performance,  emergence  of  new  facilities,  and  difficulty  In  obtaining  up-to-date 
performance  date,  such  tabulations  do  not  retain  permanent  value. 
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